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ABSTRACT 

We present a complete analysis of the Fundamental Plane (FP) of early-type galaxies (ETGs) 
in the nearby universe (z<0.1). The sample, as defined in paper I, comprises 39,993 ETGs 
located in environments covering the entire domain in local density (from field to cluster). We 
derive the FP in the grizY JHK wavebands with a detailed discussion on fitting procedure, 
bias due to selection effects and bias due to correlated errors on the effective parameters, 
and < ji >c, as key factors in obtaining meaningful FP coefficients. Studying the Kormendy 
relation (KR) we find that its slope varies from g (3.44±0.04) through K (3.80±0.02) im- 
plying that smaller size ETGs have a larger ratio of optical to NIR radii than galaxies with 
larger r^- We also examine the Faber- Jackson (FJ) relation and find that its slope is similar for 
all wavebands, within the uncertainties, with a mean value of 0.1 98 ±0.007. Writing the FP 
equation as log — a log (Jq + b < fi>e +c, we find that the "a" varies from 1.38 ± 0.02 in 
g, to 1.55 ± 0.02 in K, implying a 12% variation across the grizYJHK wavelength baseline. 
The corresponding variation of "b" is negligible (b ^ 0.316), while "c" varies by ~ 10%. 
We show that the waveband dependence of the FJ and KR results from the complex variation 
of the distribution of galaxies in the face-on projection of the FP as well as by the change of 
FP coefficients with waveband. We find that "a" and "b" become smaller for higher Sersic 
index and larger axis ratios, independent of the waveband. This suggests that these varia- 
tions are likely to be related to differences in structural and dynamical (rather than stellar 
population) properties of ETGs. It is noticeable that galaxies with bluer colours and disc-like 
isophotes have smaller "b", with the effect decreasing smoothly from g through K. Consider- 
ing a power-law relation between mass-to-light ratio and (dynamical) mass, M/ L oc M"', we 
estimate gamma from the FP coefficients in grizYJHK. The 7 decreases from 0. 224 ± 0.008 
in g, to 0.186 ± 0.009 in K band. Using the 7 values, we estimate the variation of age and 
metallicity of the stellar populations present in massive galaxies per decade in stellar mass. 
This analysis shows that in the NIR the tilt of the FP is not due to stellar population's varia- 
tion, and that ETGs have coeval stellar populations with an age variation of a few percent per 
decade in mass, and a corresponding metallicity increase of ^23%. We also show that current 
semi-analytical models of galaxy formation reproduce very well these amounts of variation 
of age and metallicity with respect to stellar mass. 
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1 INTRODUCTION 

One of the most outstanding and basic cosmological questions is 
how galaxies form and evolve. Currently the favoured scenario as- 
sumes that the assemblage of baryonic matter is driven by the evo- 
lution of dark matter haloes (Gott & Rees 1975). Given the diffi- 
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culty of observing dark matter, we rely on the luminous counter- 
part to be a beacon illuminating their evolution. The vast majority 
of stars and metals produced during the evolution of galaxies were 
formed and still reside in them. Therefore, examining the star for- 
mation history and measuring the metal content of galaxies may 
tell us how these systems evolve through cosmic time. 

The study of the global properties of elliptical galaxies took a 
substantial step forward with the application of multivariate analy- 
sis, revealing potentially meaningful scaling relations like the Fun- 
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damental Plane (FP, Brosche 1973). However, the importance of the 
technique pioneered by this paper was not immediately realised by 
the astronomical community. Determining which dimensions are 
statistically significant in a given data set is not a simple task, but 
it can reveal useful correlations involving the quantities defining a 
minimal manifold and provide insights into the physical nature of 
such correlations. Such is the case when the observed FP is associ- 
ated with the virial theorem (Djorgovski & Davis 1987; Dressier et 
al. 1987). 

Many studies over the past twenty years have tried 
to interpret the physical meaning of the FP (e.g. Faber et 
al. 19 87; Djorgovski & de Carvalho 1990; Pahre et al. 1998b; 
|j0rgensen. Franx, & Kiaer gaard 1996; J0rgens en et al.lll99 9; Dan- 
tas et al. 2003; ISernardi et aL,2003a»h.Ci; Nelan et al. 2005; Cap- 
pellari et al. 2006). The striking feature of the FP is its nar- 
rowness, implying a regularity among the global properties of 
early-type galaxies. The quantities containing the entire variance 
of the data are: effective radius, rc, central velocity dispersion, 
CTo, and mean surface brightness measured within the effective 
radius, /ic- The best representation of the FP is rc ~ f^^Io, 
where Je is mean surface brightness in flux units. Bernardi et 
al. (2003c) show a comprehensive table listing the most important 
papers presenting values of A and B and their respective errors. A 
seems to vary with the passband used in the photo metric observa- 
tion, while B does not (see e.g . Pahre, de Carvalho. & Diorgovskil 
1998b; Scodeggio et al.l [ l998l ; iMobasher et al. 1999). However, 
Bernardi et al, 2003 cL found only a marginally significant vari- 



ation of A in the SDSS optical passbands (see also Hyde & 
Bernardi 2009). La Barbera et al. (2008) (hereafter LBM08) also 
found a small difference in A when measured between r and K 
bands. 

Assuming that early-type galaxies are homologous systems in 
dynamical equilibrium and that velocity dispersion is related to the 
kinetic energy per unit mass we can write down expressions for 



mass (M) and luminosity (L), namely (M/L) 



In 



the case of a fully virialized system, A = 2 and B = -1, implying 
a constant mass-to-light ratio. However, A and B are found to dif- 
fer significantly from the virial values, resulting in the so-called 
tilt of the FP. In this case, M/L ~ M^, where gamma is ~0.25 
(Faber et al. 1987). This dependence of the mass-to-light ratio on 
galaxy mass has been interpreted as arising either from differences 
in the stellar populations or the dark matter fractions among ETGs. 
It is important to emphasize that another option to explain the tilt 
is related to the assumption that ETGs are truly virialized systems 
- in which case they should have self-similar density distributions 
and similar orbital distributions. Any departure from either or both 
of these conditions may well explain the tilt, and several studies 
have tried to disentangle these effects. For instan ce, non-homology 
seem s to contribute to at least part of the tilt jHiorth & MadsenI 
1 19951 Capelato, de Carvalho & Carlberg 1995; Ciotti, Lanzoni & 
Renzini 1996; Ciotti & Lanzoni 1997; Graham & Colless 1997; 
Busarello et al. 1997; Bertin, Ciotti, & del Principe 2002; Trujillo, 
Burkert & Bell 2004). Even studies trying to explain the tilt as a 
stellar population effect concluded that non-homology may play a 
significant role in determining the tilt of the FP (e.g. Pahre, Djor- 
govski & de Carvalho 1998b; Forbes & Ponman 1999). Another 
interesting finding from the simulations of Capelato et al. (1995) 
is that when measuring the structural parameters defining the FP 
inside larger apertures, of order a few rc, the coefficients are simi- 
lar to those implied by the virial theorem. More recently, Bolton et 
al. (2008) find a similar result when using the surface density term 
defined by the mass measured through strong lensing, and conclude 



that the tilt of the FP is due to the fraction of dark matter inside one 
effective radius (see also Tortora et al. 2009). 

This is the second paper of a series analysing the properties 
and the scaling relations of ETGs as a function of the environ- 
ment where they reside. The Spheroids Panchromatic Investiga- 
tion in Different Environmental Regions (SPIDER) utilises opti- 
cal and Near-Infrared (NIR) photometry in the grizYJHK wave- 
bands, along with spectroscopic data, taken from the UKIRT In- 
frared Deep Sky Survey-Large Area Survey (UKIDSS-LAS) and 
the Sloan Digital Sky Survey (SDSS). The selection of ETGs for 
this project is detailed in Paper I, and we refer the reader to that pa- 
per for all the details of sample selection and the procedures used 
to derive the galaxy parameters. 

In this work we focus on the derivation of the FP in the 
grizYJHK wavebands for the entire SPIDER sample. Although 
our sample contains ETGs over the entire domain of local density 
(from field to clusters), we postpone the study of the environmen- 
tal dependence of the FP to another paper in the SPIDER series 
(paper III). In the present work, we discuss the main pitfalls of the 
FP fitting procedure and how to account for selection effects and 
different sources of biases. We analyse the edge-on and face-on 
projections of the FP, as well as the two other projections of the 
FP, i.e. the Kormendy and Faber-Jackson relations. The analysis of 
these two scaling relations serves as a reference at the local Uni- 
verse (z<0.1), and at different wavebands, for other studies lacking 
data for a full FP analysis. We find a consistent picture connecting 
the waveband variation of the edge- and face-on projections of the 
FP with that of the Kormendy and Faber-Jackson relations. Finally, 
we show how the optical and NIR FPs can constrain various sce- 
narios for galaxy formation and evolution, by using the wavelength 
dependence of the FP to infer the variation of stellar population 
parameters along the ETG's sequence. 

The layout of the paper is as follows. Sec. 2 shortly describes 
the SPIDER dataset. Sec. 3 presents the different subsamples of 
ETGs used to derive the FP in grizYJHK and to analyse the 
impact of different biases on the FP. Sec. 4 details the FP fitting 
procedure. Sees. 5 and 6 analyse the Kormendy and Faber-Jackson 
relations, respectively. Sec. 7 presents one main result of this study, 
i.e. the dependence of FP slopes on waveband, from g through K. 
Sec. 8 analyses the waveband variation of the edge- and face-on 
projections of the FP. Sec. 9 describes how the optical and NIR 
scaling relations of ETGs constrain the variation of stellar popu- 
lation properties along the FP. Discussion follows in Sec. 10. A 
summary is provided in Sec. II. 

Throughout the paper, we adopt a cosmology with Ho = 
75kms-^ Mpc-\ ^2^ = 0. 3, and Qa^O-T. 



2 DATA 

The SPIDER data-set is based on a sample of 39, 993 ETGs 
(see paper I for details), with available griz photometry and 
spectroscopy from SDSS-DR6. Out of these galaxies, 5, 080 
objects have also photometry available in the YJHK wave- 
bands from UKIDSS-LAS. All galaxies have two estimates 
of the central velocity dispersion, one from SDSS-DR6 and 
an alternative measurement obt ained by fitting SDSS sp ectra 
with the software STARLIGHT jCid Femandes et al.l 120051) , us- 
ing a linear combination of simple stellar population models 
(rather than single templates as in SDSS) with different ages 
and metallicities. In both cases, STARLIGHT and SDSS-DR6, 
the (To's are aperture-corrected to an aperture of re/8, follow- 
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Table 1. Magnitude limits in grizY J H K adopted to derive tlie FP. 



waveband X" Mx Umit Nx 



g 


-19.75 


-19.71 


4467 


r 


-20.60 


-20.55 


4478 


i 


-21.02 


-20.99 


4455 


z 


-21.34 


-21.22 


4319 


Y 


-22.03 


-21.95 


4404 


J 


-22.55 


-22.54 


4317 


H 


-23.22 


-23.21 


4376 


K 


-23.60 


-23.60 


4350 



X is the equivalent magnitude limit as used in the colour-selected 
samples 

ing |j0rgensen. Franx. & Kjeergaar . In order to make 

proper comparisons to earlier studies (e.g. Bernardi et al. 2003a), 
we use SDSS velocity dispersion measurements to examine the 
scaling relations presented in this paper. In paper I, we find that the 
mean difference between (To(SDSS-DR6) and (jq (STARLIGHT) 
does not change significantly with ctq- Therefore, we do not expect 
that the choice of a given velocity dispersion measurement might 
have a dramatic impact on the FP relation. This is further discussed 
in Sees. 6, 7, and 9. 

In all wavebands, structural parameters - i.e. the effective ra- 
dius, rc, the mean surface brightness within that radius, < fi >c, 
and the Sersic index, n- have been al l homogeneously measured 
by 2DPHOT jLa Barbera et all l2008 a'). In the optical (griz), al- 
ternative estimates of the effective parameters, Vc and < fi >c, 
are also available from the SDSS-DR6 Photo pipeline. In paper 
I, we compare the different estimates of photometric and spectro- 
scopic parameters, deriving also an estimate of the 95% complete- 
ness limit of the sample in all wavebands. We find that 2DPH0T 
total magnitudes are brighter than SDSS model magnitudes, with 
the difference amounting to ~ 0.2 mag in r-band, for the faintest 
galaxies in the sample. This difference is due to the use of Ser- 
sic (2DPH0T) rather than de Vaucouleurs (Photo) models to fit the 
light distribution of ETGs, as well as to the sky estimate bias affect- 
ing SDSS effe ctive parameters jAdelman-McCarthv et alj I2OO8I : 
lAbazaiian et al . 2009). Hence, the completeness limit of the sample 
is also dependent on the source of effective parameters (2DPH0T 
vs. Photo). In r band, the sample is 95% complete at —20.32 and 
—20.55 for the SDSS and 2DPH0T parameters, respectively. In the 
following, unless explicitly said, we refer to 2DPH0T total magni- 
tudes. 



3 THE SAMPLES 

The waveband dependence of the FP is analysed using different 
subsamples of ETGs, extracted from the SPIDER sample. Details 
on each subsample are provided in Sec. 13. II In order to analyse the 
effect of different sources of bias on the FP relation, we also utilise 
several samples of ETGs, with effective parameters in r band. We 
describe the characteristics of these samples in Sec. 13.21 referring 
to them, hereafter, as the control samples of ETGs. 

3.1 The grizYJHK (SDSS-i-UKIDSS) samples of ETGs 

In order to analyse how different selection procedures might af- 
fect the dependence of the FP relation on waveband, we derive the 



FP in the grizYJHK wavebands for ETG's SPIDER subsamples 
defined by two different selection procedures. In case (i), we de- 
rive the FP for the same sample of ETGs in all wavebands, by se- 
lecting those galaxies brighter than the r-band completeness limit 
^0.07 _20.55). We exclude galaxies whose Sersic fit, in one 
of the available wavebands, has an high reduced value (^ 3). 
This cut removes less than 2% of galaxies, resulting in a sample 
of 4, 589 ETGs. In case (ii), we select different samples of ETGs 
in the different wavebands, but according to equivalent magnitude 
limits. The equivalent magnitude limits are derived by using the 
optical-NIR colour-magnitude relations (see sec. 4 of paper I). To 
this effect, we first fix the r-band magnitude limit to —20.6 and then 
translate it into the other wavebands using the colour-magnitude re- 
lations. The value of —20.6 is chosen so that, for each band, the 
equivalent magnitude limit is brighter than the completeness mag- 
nitude in that band, as defined in paper I. This makes the samples 
magnitude complete in all wavebands. The equivalent magnitude 
limits are reported in column (2) of Tab. [T] along with the 95% 
completeness magnitude limits, from paper I, in column (3), as well 
as the number of ETGs selected in each band in column (4). 

In the following, we refer to the ETG sample of case (i) as the 
(r-band) magnitude-selected sample of ETGs, while the samples of 
case (ii) are referred to as the colour-selected samples of ETGs. 

3.2 Control samples of ETGs 

We use five control samples of ETGs selected from SDSS-DR6, 
with photometry available in r band. The control samples consist 
of ETGs selected in different redshift ranges, with effective param- 
eters and central velocity dispersions measured with different meth- 
ods. In all cases, velocity dispersions are corrected to an aperture of 
re/8, following J0rgensen. Franx. & Kiaergaard Each con- 

trol sample is named with a letter, as shown in Tab. (2] where we 
summarise the basic characteristics of the five samples. 

- Sample A is obtained from the sample of 39, 993 ETGs defined 
in paper I. We select all galaxies with an r-band model magnitude 
brighter than —20.32. This magnitude cut corresponds to the 95% 
completeness limit in r band, as defined in paper I, when using 
SDSS model magnitude s (see Sec.[2l(. Effective parameters are ob- 
tained from SDSS, as in lBemardiet^j2003al) . 

- Sample B is a subsample of sample A, consisting of all the 
ETGs that also have photometry available in the YJHK wave- 
bands (see paper I). Such sample is used to estimate the impact of 
matching SDSS to UKIDSS data on the FP relation. 

- Sample C is defined to explore a wider magnitude range than 
that of samples A and B. We query the SDSS-DR6 database for 
ETGs in a redshift range of 2: = 0.02 to 0.03. ETGs are defined 
according to the same criteria as in paper I, i.e. zwarning = 0, 
eclass < and fracDeVr > 0.8. No requirement is done for 
the galaxy velocity dispersion. The query results into a list of 3732 
galaxies, that hereafter we refer to as the low-redshift sample of 
ETGs. All galaxies have effective parameters from SDSS. Using 
the same procedure as in paper I, we estimate a 95% completeness 
limit of ""^Mr = —17.64 (model magnitude). Since velocity dis- 
persions from SDSS are not available for all galaxies in this sample, 
we assign fake ao values to each pair of r^ and < >e values, as 
described in Sec. 14.21 

- Samples D and E are defined in the same way as samples A 
and B, respectively, but using 2DPHOT rather than SDSS effective 
parameters. For both samples, we select all the ETGs in the SPI- 
DER sample with total magnitude brighter than ""^A'/r= —20.55 
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(corresponding to the 2DPHOT completeness magnitude). Sam- 
ple E is obtained from sample D by selecting only those objects 
with matched photometry in UKIDSS. Sample E coincides with 
the magnitude-selected sample of ETGs in r band (Sec. l3.lt . 



4 DERIVING THE FP 

We write the FP relation as: 

logrc = alogCTo + b <M>c +c, (1) 

where "a" and "b" are the slopes, and "c" is the offset. We de- 
note the rms of residuals around the FP with respect to log as 
Sre, referring to "a", "b" , "c", and s,-^ as the coefficients of the 
FP. We estimate the FP coefficients by a procedure consisting of 
three steps. First, we derive the values of "a" and "6", as described 
in Sec. 14.11 The slopes are then corrected for different sources of 
biases, including selection effects (Sec. 14.21 ) and the effect of cor- 
related uncertainties on log re and < /i >e, (Sec. I4.3l >. The bias- 
corrected values of "a" and "6" are then used to estimate "c" and 
(Sec. 14. id . This procedure is tested through the ETG's control 
samples, as discussed in Sec. 14.41 

4.1 Fitting procedure 

We obtain a first estimate of "a" and "b" by minimising the sum 
of absolute residuals around the FP. When compared to the ordi- 
nary least-squares fitting method, where one minimises the sum 
of squared residuals, this procedure is more robust, being less 
sensiti ve to outliers in the distribution of da ta-points around the 
plane ('j0rHensen, Franx, & Kisreaar 3 '199^ (hereafter JFK96). 
We adopt two different fitting methods, by minimising the resid- 
uals in log (Jo and the orthogonal residuals about the plane. The 
orthogonal fit - adopted in most of previous works - has the main 
advantages of treating all the variables symmetrically, while the 
log (To regression is essentially independent of selections effects 
in th e plane of effective parameters, such as the magnitude limit 
(see iLa Barbera. Busarello. Capacciolil [lOOOl . hereafter LBCOO). 
The values of "a" and "b" are corrected for selection effects and 
correlated errors on effective parameters (see Sees. 14.21 and 14.3b . 
The value of "c" is then derived as the median value of the quan- 
tity logr^ — a log (To — b < fi >e, over all the galaxies of a given 
sample, with "a" and "b" being the bias-corrected values of FP 
slopes. As shown in Sec. 4.2, when compared to the more common 
practise of estimating "c" through the least-squares procedure it- 
self, the above estimate has the advantage of providing an unbiased 
value of "c" , regardless of the magnitude selection of the sample. 
For both fitting methods, we calculate the scatter of the FP, Sr^, 
from the mean value of the absolute residuals in log around the 
plane, using the bias-corrected slopes. As for "c", this procedure 
provides an unbiased estimate of the FP scatter (see Sec. |4.2t . 

4.2 Bias due to selection effects 

To estimate how selection criteria (e.g. the magnitude limit) affect 
the FP coefficients, we use a simulated sample of data-points in the 
space of logre, < fi >c, and logtTo, resembling the distribution 
of ETGs in that space. The simulated sample is created from the 
control sample C, namely all ETGs from SDSS-DR6 in the redshift 
range of 0.02 to 0.03, brighter than an r-band model magnitude of 
""^Mr = —17.64 (Sec. |3.2t . Since galaxies in this sample do not 
have available velocity dispersions, we assign fake ao values. For 



each galaxy, we use its log and < /i >e to obtain a value of ao 
from the FP relation (Eq. [TJ. That value is then shifted according 
to a random Gaussian deviate, with a given width value So-q, that 
describes the scatter of the FP along the ao axis. The slopes, offset, 
and scatter parameters are chosen with an iterative procedure. 

- First, we select all galaxies in sample C with available ao from 
SDSS-DR6, applying similar cuts in magnitude and velocity dis- 
persion as those for the r-band magnitude-selected sample of ETGs 
(sec Sec. 13.1b . This is done by selecting all galaxies with model 
magnitude brighter than -20.280 and 70 ^ (tq ^ 420 fcms"^ 
This subsample consists of 1682 ETGs out of 3690 galaxies in sam- 
ple C. We derive the best-fitting FP coefficients for this subsample, 
referring to them as the reference coefficients of the FP. 

- We assign fake ao values to sample C by using guess values of 
"a", "6", "c" and Sc„. Applying the same cuts in magnitude and 
velocity dispersion as in the above step, we derive the best-fitting 
FP coefficients and compared them to the reference FP coefficients. 
The guess values of "a", "6", "c" and So-q are changed until the 
best-fitting simulated FP matches the reference relation. In practise, 
we are able to match the simulated and reference coefficients at 
better than 2% for both the log ao and orthogonal regressions. 

Fig. [T] compares the distribution of the 1682 ETGs with available 
(To's from SDSS in sample C with that of data-points for one of the 
toy samples, showing the similarity of the two distributions. The 
above procedure allows us to create simulated samples in the space 
of logTc, < /i >c, and log (To down to a (model) magnitude limit 
of —17.64, which is more than 2.5 magnitudes fainter than the r- 
band completeness limit (—20.32) of the ETG samples of Sec. 13. II 
The effect of any selection cut on the FP can then be estimated by 
computing the relative variation of FP coefficients as one applies 
that selection to the toy samples. 

Fig.[2]plots the relative variation of FP coefficients as a func- 
tion of the magnitude cut. The relative variation of a given quantity, 
X, out of "a", "6", "c", and s^,, is computed as [xcut — x)/x, 
where Xcut is the value estimated for that quantity when the cut is 
applied. Here, instead of using the procedure of Sec. 14. 11 the value 
of "c" is directly derived from the fit, and the s,-^ is obtained as 
the mean absolute deviation of residuals around the plane, using 
the (no bias-corrected) best-fitting coefficients "a", "fe", and "c". 
For the orthogonal fit, we see that brighter the magnitude cut, more 
the FP coefficients tend to be underestimated. T his finding is con- 
sisten t with that of previous studies (see LBCOO. lHvde & Bernardil 
'200^. For the log ao fit, the FP coefficients are very insensitive, 
as somewhat expected, to the selection in magnitude. The vertical 
lines in Fig. [2] correspond to the (r-band) model magnitude limit 
(^■°''Mr^ -20.28) of the magnitude- selected sample of ETGs 
(see Sec. I3.lt , after the small amount of luminosity evolution be- 
tween z = 0.025 and z = 0.075 has been removed (see above). 
For that magnitude limit, the amounts of bias in "a", "c" 
and (horizontal lines) are significant, amounting to about 27%, 
8%, 16% and 22%, respectively. The same amounts of bias are 
also expected to affect the colour-selected samples of ETGs (see 

^ Notice that the value of —20.28 is 0.04 mag fainter than the r-band 
model magnitude limit of the magnitude-selected ETG sample (^■'^'^ Mr = 
-20.32, see Sec.|2) We obtain -20.28 by adding to -20.32 the difference 
of evolutionary correction between the median redshift of the ETG's sam- 
ple of pap er \(z = 0.0725) and t hat of sample C (z = 0.025). To this aim, 
following 'B emardi et al (hereafter BER03b), we parametrize the 

evolutionary correction as — 2.5Qlog(l -I- z), where the coefficient Q is 
equal to ~ 0.85 in r-band , at redshift z < 0.3. 
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Table 2. Control samples of ETGs. 





A 


B 


C 


D 


E 


Number of galaxies 


37273 


4796 


3690 


36205 


4589 


Redshift range 


0.05 ^ ^ ^ 0.095 


0.05 ^ ^ ^ 0.095 


0.02 ^ 2 s£ 0.03 


0.05 ^ 2 s£ 0.095 


0.05 ^ 2 ^ 0.095 


Limiting ""^Mr 


-20.32 


-20.32 


-17.64 


-20.55 


-20.55 


Source of ro and < ^ >c 


SDSS 


SDSS 


SDSS 


2DPHOT 


2DPHOT 


Available wavebands 


r 


grizYJHK 


r 


r 


grizYJHK 



CD 
tN 



o 
b 

Oi CM 



CM - 



a=1.56±0.05 
b=0.341 ±0.006 



=0.067±0.002 



a=1.54±0.05 " 
b=0.344±0.007 . 

=0.068±0.002 
I 1 I c 



o 



3.2 



3.4 



3.6 



3.8 




CO 




19 -18 



log r - 0.34 </J.>- + 9.7 



0.07 



0.07 



Figure 1. Comparison of the short edge-on projection of the FP for galax- 
ies in sample C with available erg from SDSS (black) and one toy sam- 
ple (red). For both samples, only points with ^■^'^ ^5 —20.28 and 
70 ^ (TO ^ 420 km/s have been selected. The values of log ctq are plotted 
against the variable logre— 6< /j. >c, i e. the combination of photometric 
parameters entering the FP. The slopes ("a" and "b") and scatter (Sr,) of 
the relation, obtained from the log ctq regression procedure, are reported in 
the upper-left and lower-light corners of the plot for the observed and toy 
samples, respectively. Notice the negligible difference between the two sets 
of coefficients. 



Figure 2. Relative variation of FP coefficients as a function of the mag- 
nitude cut (see the text). The variation is computed between the magni- 
tude selected and entire toy samples. Empty and filled circles con'espond 
to the results obtained for the log ag and orthogonal fits, respectively, as 
shown in lower-right comer of the upper-left panel. From left to right and 
top to bottom, the four panels show the relation variation (bias) in "a", 
"b", "c", and Sr^, respectively. The vertical and horizontal dashed lines 
mai'k the completeness of the magnitude-selected sample of ETGs and the 
corresponding bias values, respectively. The red dashed curve in each panel 
is the fourth order polynomial fit performed to model the bias as a function 
ofOO^M,. 



Sec. 13. lb . whose magnitude limit in r-band is very similar to that of 
the magnitude-selected sample. We also used the simulated sam- 
ples to estimate the impact of the ctq cut of the ETG's sample on 
the FP relation. To this aim, we selected only simulated points with 
magnitudes brighter than ° °^Mr= —20.28. Applying the ctq se- 
lection (70 ^ (To ^ 420 km/s), we found that relative variation of 
FP slopes is completely negligible (<1%). This is due to the fact 
that, for the magnitude range considered here, almost all galaxies 
have (To > 70 km/s, making the ctq selection unimportant. 

For each sample of ETGs, as defined in Sec. 13.11 we con- 
sider the corresponding 2DPHOT r-band magnitude limit. For the 
magnitude- and colour-selected subsamples, these limits amount to 
—20.55 and —20.6, respectively. The 2DPH0T magnitude limit 
is translated to a (model) magnitude limit by adding the term 
0.23 mag which is the difference of 2DPHOT and SDSS complete- 
ness magnitudes (Sec.|2l(. For a given sample, the amount of bias 
on "a" and "b" is then estimated evaluating the trends in Fig.|2]for 
the r-band model magnitude limit of that sample. This is done only 
for the orthogonal regression procedure, by modelling the trends in 



Fig.|2]with fourth order polynomials. The biased values of "a" and 
"b" are multiplied by the estimated x/xcut factors. Notice that the 
same correction factor is applied to all the grizYJHK wavebands 
(see also Sec. l7.lt . The bias-corrected values of "a" and "b" are 
used to estimate "c" and Sr^ (Sec. l4.lt . Fig.[3]shows how the val- 
ues of "c" and Sr^ vary as a function of the magnitude limit, when 
this procedure is applied, rather than estimating c and Sr^ from the 
fit, as in Fig.|2] As stated in Sec. 14.11 the estimates of "c" and s,.^ 
from the bias-corrected values of "a" and "b" are almost insensi- 
tive, within ~ 2%, to the magnitude selection. 

4.3 Bias due to correlated errors on re and <ii>c 

Another possible source of bias on FP coefficients is the correla- 
tion of uncertainties on logre and < >c. As shown in paper I, 
the errors on effective parameters mainly depend on the signal-to- 
noise per pixel of galaxy images, and are slightly larger in the NIR 
than in the optical wavebands. For instance, the median value of the 
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Figure 3. Relative variation of the FP offset (upper panel) and scatter (lower 
panel) as a function of the magnitude cut, as estimated from the toy FP 
samples. The offset, "c", and scatter, s,.^, are estimated from the bias- 
con'ected FP slopes. Empty and filled circles correspond to the log ctq and 
orthogonal fit values, respectively. The red dashed lines show a fourth order 
polynomial fit of the filled circles. The vertical and dashed lines mark the 
completeness of the magnitude-selected ETG sample and the corresponding 
expected bias values. Notice that the bias is negligible for both quantities, 
being smaller than ~ 2%. 



log uncertainties increases from ~ 0.09 in g-band, to ~ 0.14 in 
K-band. This variation might imply a spurious dependence of FP 
coefficients on waveband, and thus we have to correct the FP slopes 
separately in each band. 

The corrections are estimated by (1) constructing simulated sam- 
ples of data-points in the space of logrc, < fi >c, and logao, 
resembling the distribution that galaxy's parameters would have 
in that space if no correlated errors on and < /i >c would be 
present (Sec. |4.3.Tt . and (2) estimating how the FP slopes change 
by adding correlated uncertainties on the effective parameters of 
such simulated samples (Sec. I4.3.2t . Notice that the toy samples 
of Sec. 14.21 are not suitable to apply the above procedure, since 
the corresponding effective parameters already include the effect 
of correlated errors on the effective parameters. 



4.3.1 Simulated samples with no correlated errors 

Each simulated sample is generated as follows. First, we extract 
log Te values from a random deviate whose centre and width values 
are given by the mean (0.27 dex) and standard deviation (Se = 
0.25 dex) of the log distribution of sample C. For a given log r^, 
we assign a < >o value by the Kormendy relation (hereafter KR) 



where p\ and p2 are the offset and slope, respectively. The values 
of pi and p2 are derived by a robust least-squares fitting procedure 
for galaxies in sample C, by minimising the absolute sum of resid- 
uals in < >c around the relation. As shown bv La Barbera et al.l 
{2003) (hereafter LBM03), the KR fit is quite insensitive to the cor- 
related errors on log and <p,>e. The fit gives pi ~ 18.969 and 
P2 ~ 1.95, respectively. Then, we shift the values of < >c ac- 
cording to a normal Gaussian deviate of width 0.4:mag/arcsec^, 
corresponding to the intrinsic dispersion in < ^ >c of the KR 
(LBM03). For a given pair of logrc and < ^ >c values, we as- 
sign a log (Jo value by the FP relation (Eq. [T]!. The logao values 
are shifted according to a random deviate with given width, sq. The 
free parameters of this procedure, i.e. the FP slopes and offset, and 
the value of so, are chosen so that, on average, the FP coefficients 
of simulated samples match those of the magnitude- selected sam- 
ple of ETGs, with the same iterative procedure as in Sec. 14.21 



<fI>e=Pl +P2 log re 



(2) 



4.3.2 The effect of correlated uncertainties 

The logro and < fj,>c of the simulated samples are then shifted 
according to a two-dimensional random deviate, whose covariance 
matrix terms are given by the median uncertainties on log rc and 
< /I >o for galaxies in the magnitude-selected samples of ETGs. 
The procedure is repeated for each waveband, by using the cor- 
responding median covariance matrix of uncertainties on effective 
parameters. We derive the FP slopes by (i) applying the correlated 
errors, and (ii) without applying any simulated uncertainty on the 
effective parameters. We indicate as Sa and St the ratios of FP 
slopes of case (ii) with respect to those obtained in case (i). Each toy 
sample includes A'^ = 2000 data-points, and the values of 5a and 
5t are averaged over 300 realisations. The values of 5a and 5t in 
grizY J H K bands, for both the orthogonal and log ao regression 
procedures, are reported in Tab|3] The correlated uncertainties on 
effective parameters tend to increase the value of the log ao slope 
of the FP, and decrease the coefficient of the < ^ >c term. The 
effect is quite small, in particular for the coefficient "a" , amount- 
ing to less than a few percent. The bias is larger for "fe", and varies 
almost by a factor of two from the optical to the NIR wavebands. 
Moreover, unlike the bias due to selection effects, it affects both 
the orthogonal and log ctq regression procedures. Due to the large 
number of galaxies in the SPIDER sample, the factors in Tab.[3]are 
not negligible with respect to the typical errors on FP slopes (see 
Sec. 17. It . Hence, we correct the slopes of the FP in each band mul- 
tiplying them by the corresponding 5a and 5b factors in Tab [3] We 
have also performed some tests to check how robust the values of 5a 
and 5i are with respect to the procedure outlined above. First, one 
can notice that the adopted slope of the KR (p2 ~ 1.95) is smaller 
than that of p2 ~ 3 found by other studies (see LBM03 and ref- 
erences therein) and by that reported for the SPIDER samples in 
Sec. |5] Hence, we derived the offset of the KR by fixing p2 = 3 
and repeated the above procedure with the corresponding values of 
pi and p2- Second, one may notice that the width value itself of the 
log rc distribution, Se = 0.25 dex (see above), is broadened by the 
measurement errors on logrc, and hence does not correspond to 
the intrinsic width of the log ro distribution. To account for this ef- 
fect, we repeated the above procedure by subtracting in quadrature 
0.1 dex (the typical uncertainty on log rc in r-band) to the value of 
Se. For both tests, we found that the variation of the 5a and 5b esti- 
mates in Tab.[3]is completely negligible, being smaller than 0.5%. 
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Table 3. Effect of the correlated uncertainties of effective parameters on 
tlie slopes of the FP in different wavebands. 



orthogonal fit log ctq fit 
waveband 5a &h &h 



9 


0.995 


1.038 


0.984 


1.039 


r 


0.990 


1.035 


0.980 


1.035 


i 


0.996 


1.031 


0.992 


1.031 


z 


0.986 


1.043 


0.975 


1.043 


Y 


0.992 


1.040 


0.980 


1.039 


J 


0.985 


1.060 


0.977 


1.062 


H 


0.980 


1.066 


0.964 


1.065 


K 


0.975 


1.070 


0.956 


1.069 



4.4 Comparison of bias-corrected FP coefficients in r-band 

In order to test the above procedure for deriving the coefficients of 
the FP and correct them for the different sources of biases, we ap- 
ply it to the control samples of ETGs (Sec. l3.2t . In Fig. ID we plot 
the corrected slopes of the r-band FP for the five control samples. 
For sample C, we select only those 1682 galaxies with available 
ao's from SDSS, and (model) magnitudes brighter than —20.28 
(see Sec. 14. 3t . The values of "a" and "h" are also c ompared to 
those recently obtained from i Hvde & Bemardil ( |2009|) . who took 
into account selection effects in the fitting procedure, rather than 
applying correction factors as we do here. For all samples, the FP 
slopes are corrected for the magnitude bias evaluating the polyno- 
mial curves in Fig.|2]at a (model) magnitude of Mr= -20.28. 
For samples D and E, this ""^Mr value corresponds to the magni- 
tude limit of —20.55, after difference between model and 2DPH0T 
total magnitudes is taken into account (see paper I). In order to re- 
move the effect of correlated errors on effective parameters, the 
slopes of samples D and E have also been divided by the r-band 
correction factors reported in Tab. [3pl 

Fig. |4] shows that the FP slopes of all control samples are re- 
markably consistent wit hin the 2a level, and differ by less than 
~ 3% from the values of iHvde & Bernardil ( [2009i) . proving the ro- 
bustness of the procedure outlined above to derive bias-corrected 
FP coefficients. The consistency of FP slopes between samples A 
and B (D and E) shows that matching the ETG's sample with 
UKIDSS does not lead to any significant bias in the estimate of 
FP coefficients, in agreement with LBM08. One can also notice 
that, although the SDSS and 2DPHOT effective parameters differ 
significantly (see paper I), the corresponding FP relations are very 
consistent, as shown by the consistency of FP slopes between sam- 
ple A and D (B and E). This is due to the fact that the combination 
of To and < >o that enters the FP is determined with much better 
accuracy that Vc and < ^ >c themselves (see .Kelson et al..,200Q) . 
making the FP relation very stable. 



5 THE KORMENDY RELATION 

Fig.|5]plots the rc-< fi >c diagram for the colour- selected sam- 
ples of ETGs (Sec. I3.lt . from g through K. For each band, the 

^ We also estimated the uncertainties on SDSS Photo parameters in the 
same way as for the 2DPH0T effective parameters, i.e. by comparing the 
values of re and < >e from SDSS in r and i bands (see paper I). For these 
uncertainties, we found that the r-band correction factors on FP slopes are 
even smaller than those reported in Tab. [3] Hence, we decided not to apply 
any further correction factor to samples A, B, and C. 
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Figure 4. Slopes of the r-band FP, corrected for selection effects and cor- 
related errors on effective parameters, for the control samples of ETGs 
(Tab. |2). Each sample is plotted with a different colour, as shown in the 
upper-left corner For each point, the corresponding concentric ellipses de- 
note the one and two a confidence contours for a two-dimensional normal 
Gaussian d eviate. T he dashed lines mark the values of "a" and "6" ob- 
tained from lHvde & Bemardil |2009|) . 



figure also exhibits the completeness limit of the sample in that 
band, from Tab. [T] Galaxies follow a well-defined KR in all wave- 
bands. We write the KR as in Eq. |2] In order to characterise the 
offset, pi, the slope, p2, and the scatter, s^,^, of the KR, we ap- 
ply the modified least-squares (hereafter MLS) fitting procedure 
of LBM03. The MLS fit allows the coefficients of the KR to be 
derived by accounting for selection cuts in the r^-K fj.>e diagram, 
such as the magnitude limit. LBM03 applied three MLS fits. The 
MLSiogr^ and MLS<jj>^ regressions are obtained by minimising 
the residuals around the relation with respect to log re and < /i >c, 
respectively. The MLSB fit corresponds to the bisector line of the 
MLSiog,-^ and MLS<fi>^ fits. The MLSB method is more robust 
and effective (i.e. lower uncertainties on fitting coefficients) with 
respect to the other MLS fits. For this reason, we apply here only 
the MLSB fit. Moreover, we generalise the MLS method to the case 
where orthogonal residuals around the relation are minimised. This 
orthogonal MLS fit (hereafter MLSO) is described in App. |A] For 
both the MLSB and MLSO fits, the KR coefficients are derived ac- 
counting for the magnitude limit of the sample in the corresponding 
waveband. The scatter of the KR is obtained by the standard devia- 
tion of the log r-c residuals about the line, accounting for the mag- 
nitude cut as detailed in App. IaI Fig. [5] also plots the MLSB and 
MLSO lines. The corresponding fitting coefficients are reported in 
Tab.|4] 

From Tab.|4]one sees that the MLSO fit gives a larger value of 
the slope, p2, with respect to the bisector fit. The scatter around the 
KR is independent of waveband, and larger, by ~ 0.01 dex, for the 
MLSB than for the MLSO fit. The KR smoothly steepens from the 
g through the K band. This is shown in Fig. |6] where we plot the 
MLSB slope of the KR as a function of the logarithmic effective 
wavelength of each filter The p2 smoothly increase from a value of 
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Figure 5. Kormendy relation of ETGs in tlie grizY J H K wavebands (from left to right and top to bottom). For each panel, the dotted line mark the 
completeness magnitude in the corresponding waveband. The solid and dashed lines are the best-fitting relations, obtained by the orthogonal and bisector 
fitting methods, respectively (see the text). 
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Figure 6. The slope of the KR, obtained with the MLSB fit, p2 , is plotted as 
a function of the logarithmic effective wavelength, log A, of the passbands 
where effective parameters are measured. 



~ 3.44 in g to 3.8 in K. A similar trend is also observed for the 
results of the MLSO fit. 

In order to analyse the trend of p2 with wa veband, we follow 
the same approach as in iLa Barbera et^ j2004h . Given two wave- 
bands, X and W (X, W = grizYJHK), one can relate the corre- 
sponding slopes of the KR, P2,x and P2,w, through the following 
equation: 



{P2,w + 5Axw) — C(5 — P2,w) 
1 + Axw 



(3) 



where Axw is the slope of the X — W ws. W colour-magnitude 
relation, and ( parametrizes the variation of the mean logarithmic 
ratio of X to W effective radii, los , as a function of Vc'. 



, re,w ^, 

log oc Clogre,x- 



(4) 



We first consider the case where the mean ratio of effective radii 
does not change along the sequence of ETGs, i.e. ^ = 0. Set- 
ting X — K and W = g, using the value of AgK from paper I 
(0.034 ± 0.016) and the value of p2,g of the MLSB fit, Eq.[3]would 
imply P2,K = 3.53 ± 0.02. This value is significantly smaller than 
that reported in Tab.|4](3.80 ± 0.02), implying that the assumption 
^ = is incorrect. Indeed, inverting Eq. [3] and using the MLSB 
values of p2.g and P2,K, one obtains — —0.19 ± 0.02. The nega- 
tive sign of (" implies that galaxies with smaller re,K tend to have, 
on average, also larger value. In other terms, the NIR light 
profile of ETGs is more concentrated in the centre with respect to 
the optical for small (relative to larger) galaxies. The dependence 
of on re, if can be directly analysed by binning the SPIDER 
sample with respect to r^^K and computing the median value of 
in each bin. The result of this procedure is shown in Fig. [7] 

We clearly see that the median value of -^^^^ decreases as re k in- 
creases, and that the trend is fully consistent with what expected 
from the waveband variation of the KR slope (Eq.|3]and Eq.|4j see 
dashed line in the Figure). In the simplistic assumption that ^^^ji 
is a good proxy for the internal colour gradient in ETGs, the in- 
creasing of KR slope from g through K would imply that smaller 
size ETGs have stronger (more negative) internal colour gradients 
than galaxies with larger r^. This point will be further analysed in 
a forthcoming paper, studyin g the optical-NIR colour gradien ts in 
the SPIDER sample (see also iLa Barbera & de Carvalhol2009l) and 
their correlations with galaxy properties. 
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log r^^^ (kpc) 

Figure 7. Logaiithmic ratio of g to /f-band eifective radii as a function 
of log re in K-band. The solid line connects the data-points obtained by 
median binning the distribution of logre,g/re_/f with respect to log r^^K, 
with each bin including the same number (N=200) of points. The trend is 
fully consistent with that expected from the variation of KR slope from g to 
K (dashed Hne). The dotted line marks the value of zero. En'or bars denote 
la errors on median values in different bins. 



The slope of the MLSB fit can be compared to that obtained 
from LBM03 for a sample of ETGs in clusters at intermediate red- 
shifts, from 2 ~ to z ~ 0.64. Using the MLSB fit, LBM03 found 
P2 — 2.92 ± 0.08 in V-band restframe. This should be compared 
with the value of p2 = 3.44 ± 0.04 we obtain for the SPIDER 
sample in the g band (see Tab. which matches approximately 
V-band restframe. The slope of LBM03 is significantly flatter, by 
~ 15%, that that we find here. One should notice that LBM03 se- 
lected ETGs by a cut in the Sersic index n (n > 2), while ETGs 
are defined here according to several photometric and spectroscopic 
criteria. Moreover, ETGs in the SPIDER sample reside in a wide 
range of environments, while ETGs in LBM03 mostly belong to 
rich galaxy clusters. Both these issues might be responsible for the 
above difference of KR slope values. 



6 THE FABER-JACKSON RELATION 

We write the Faber-Jackson (hereafter FJ) relation as: 

logcro = Ao + Ai(logL + 0.4X) (5) 

where Aq and Ai are the offset and slope of the relation, and X 
is the magnitude limit in a given waveband (see Tab. [TJ. Accord- 
ing to this notation, the coefficient Aq is the logcro value pre- 
dicted from the FJ relation for galaxies of magnitude X. The galaxy 
luminosity, L, is defined as iQ-O-'ix" "'' m ^ ^^^^^ o.m j^,^ 
2DPHOT absolute magnitude in the given band. In order to de- 
rive the coefficients Aq and Ai we use the colour-selected samples 
of ETGs (Sec.lHB- Fig. [8] plots the distributions of ETGs in the 
logao vs. logL diagrams. Each sample is binned in logcro, and 
the peak value of the log L distribution in a given bin is co mputed 
by the bi- weight statistics l lBeers. Flvnn. & Gebhardtlll990h . Since 



all colour-selected samples are magnitude complete, the binning 
procedure produces unbiased estimates of the average log L value 
as a function of log cro. The binned values of log L vs. log ctq are 
then fitted with an orthogonal least-squares fitting procedure. For 
each band, the fit is performed over a fixed luminosity range of one 
decade, with — 0.4X ^ logL ^ — 0.4X + 1. Uncertainties on 
Ao and Ai are estimated by TV = 500 bootstrap iterations, shifting 
each time the log L binned values according to their error bars. The 
values of Ao and Ai in grizYJHK bands are reported in Tab. [6] 
along with the log cro scatter of the relation, j , and its intrinsic 
dispersion, o"^^ . The scatter is estimated as follows. For each boot- 
strap iteration, we calculate the rms of the log ao residuals through 
the median absolute deviation estimator. The mean value and the 
standard deviation of the rms values among the different iterations 
provide the j , and its error. The intrinsic scatter is computed by a 
similar procedure, subtracting in quadrature, for each iteration, the 
amount of dispersion due to the uncertainties on log L and log cro 
from the rms values. Considering the uncertainties, the slopes of 
the FJ relations are consistent among the different wavebands, 
with the mean value of Ai amounting to 0.198 ± 0.007. Using 
the magnitude- rather than the colour-selected samples of ETGs, 
this result does not change, with the value of Ai varying from 
0.192±0.018 to 0.209±0.018 in r band, and from 0.220±0.023 to 
0.216 ±0.032 in K band. Using STARLIGHT log ao values would 
also not change significantly the Ai values, with the mean value of 
Ai varying from 0.198 ± 0.007 to 0.187 ± 0.007. For what con- 
cerns the intrinsic dispersion around the FJ relation, it smoothly 
decreases by ~ 0.008 dex from g through K, with a value of 
~ 0.091 dex in the optical and ~ 0.083 dex in K band. Fixing 
the slope of the FJ relation in all wavebands to the average value of 
Ai= 0.198, would make this amount of variation to be 0.017 dex, 
rather than 0.008 dex. Subtracting in quadrature the values of ^ 
between the g- and K-bands, one obtains a value of ~ 0.037 dex 
(i.e. ~ 9% in ao). 

The slope value of the r-band FJ relati on is close, but flatter, 
than that of 0.25 reported bv lBemardil ( 12007 ") (see their eq. 2). This 
difference can be explained by the fact that we use Sersic (rather 
than de Vaucouleurs) total magnitudes and by the small systematic 
effect in SDSS model magnitudes (see paper I). As shown in pa- 
per 1, both effects make 2DPH0T total magnitudes to be shifted 
toward brighter values with respect to SDSS model magnitudes. 
The amount of shift is larger for bright than faint galaxies, pro- 
ducing a flatter FJ relation. The difference might also be related 
to the fact that the slope of the FJ relation seems to change ac- 
cor ding to the magnitude rang e where galaxies are selected (see 
e.g. lMatkovic & Guzmanll2005r) . The slope value of the K-band re- 
lat ion, Ai^ 0.23, is fully consistent w i th the v alue of 0.24 reported 
by IPahre, de Carvalho, & Diorgovskil ( Il998ah . For what concerns 
the intrinsic dispersion, we find a value of ct^ ^ ~ 0.09 dex in the 
optical, while iBernardi. t2007i) . find a smaller value of ~ 0.07 dex. 



7 FP SLOPES 

7.1 Variation from g through K 

Due to the large sample size and the wide wavelength baseline 
provided by SDSS-l-UKIDSS, we can establish the waveband de- 
pendence of the FP with unprecedented accuracy. Fig. |9] plots 
the slopes of the FP in different wavebands, obtained for the 
magnitude- and colour-selected subsamples of ETGs (Sec. 13. it . In 
each case, we show the results of both the log ao and orthogonal 
regression procedures. 




Figure 8. The Faber- Jackson relation of ETGs in the grizY J H K wavebands (from left to right and top to bottom). For each panel, the grey curve is obtained 
by binning the data with respect to log L, with each including the same number {N = 40) of galaxies. For each bin, the bi-weight peak of the log ctq 
distribution is computed. Coloured lines show the orthogonal fit to the binned data. 




b b 

Figure 9. The log ctq slope, "a", of the FP is plotted against the < fj,>c slope "b" . The left panel shows the case where the same sample of ETGs is used 
to derive the FP in all wavebands, while the right panel exhibits the results obtained for the colour-selected samples (see Sec. 13. U . In each panel, different 
colours denote different wavebands as shown in the lower-left corner of the figure. Filled and empty symbols mark the results of the orthogonal and log ctq 
fits, respectively, with dashed and solid ellipses corresponding to la confidence contours. 
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Table 4. Coefficients of the Kormendy relation in grizYJHK. 



band 


Pi 


P2 




Pi 


P2 








MLSB fit 






MLSO fit 




g 


19.16 ±0.04 


3.44 ± 0.04 


0.126 ± 0.002 


18.92 ±0.02 


3.68 ± 0.02 


0.115 ± 0.001 


r 


18.16 ±0.02 


3.55 ± 0.02 


0.120 ±0.002 


18.02 ±0.02 


3.72 ± 0.03 


0.114 ±0.001 


i 


17.74 ±0.02 


3.60 ±0.02 


0.122 ± 0.002 


17.60 ±0.02 


3.74 ±0.02 


0.117 ±0.002 


z 


17.42 ±0.02 


3.61 ± 0.02 


0.121 ± 0.002 


17.29 ±0.02 


3.73 ± 0.03 


0.116 ±0.002 


Y 


16.59 ±0.02 


3.67 ±0.02 


0.125 ± 0.002 


16.39 ±0.02 


3.90 ±0.03 


0.117 ±0.001 


J 


16.03 ±0.02 


3.72 ± 0.02 


0.126 ± 0.002 


15.84 ±0.02 


3.95 ± 0.02 


0.117 ±0.002 


H 


15.31 ±0.02 


3.77 ±0.02 


0.126 ±0.002 


15.13 ±0.02 


3.99 ± 0.03 


0.119 ± 0.001 


K 


14.91 ±0.02 


3.80 ±0.02 


0.128 ±0.002 


14.70 ±0.02 


4.04 ±0.03 


0.119 ± 0.001 



Table 5. Coefficients of the Faber- Jackson relation in grizY J H K . 



band Aq Ai o-p, a'^ 



9 


2.158 ±0.008 


0.172 ±0.018 


0.097 ±0.002 


0.091 ±0.002 


r 


2.151 ±0.008 


0.192 ±0.018 


0.096 ± 0.002 


0.090 ± 0.002 


i 


2.155 ±0.008 


0.185 ±0.016 


0.093 ± 0.002 


0.087 ± 0.002 


z 


2.158 ±0.008 


0.172 ± 0.018 


0.097 ±0.002 


0.091 ± 0.002 


Y 


2.144 ±0.007 


0.217 ±0.016 


0.094 ± 0.002 


0.087 ± 0.009 


J 


2.168 ±0.008 


0.194 ±0.022 


0.091 ±0.002 


0.084 ± 0.002 


H 


2.140 ±0.008 


0.233 ± 0.018 


0.091 ± 0.002 


0.084 ± 0.002 


K 


2.143 ±0.009 


0.220 ±0.023 


0.090 ± 0.002 


0.083 ± 0.002 



The slopes of the orthogonal fit are corrected for the magni- 
tude cut bias as described in Sec. 14.21 In the r band, the 2DPH0T 
magnitude limit of the magnitude- and colour-selected samples of 
ETGs are —20.55 and —20.60, respectively. We convert these val- 
ues to model magnitude limits in r-band at redshift z ~ 0.025, 
and then estimate the corresponding correction factors on "a" and 
"h" from the polynomial curves in Fig. |2] Since we have selected 
either the same sample of ETGs at all wavebands, or ETG's sam- 
ples with equivalent magnitude limits (i.e. colour-selected sam- 
ples, see Sec. 13. It , we apply the same correction factors to all the 
grizYJHK wavebands. Therefore, although the values of FP co- 
efficients in a given band depend on the correction factors, their 
relative variation from g through K is essentially independent of 
them. For the log ao fitting method, which is not affected from the 
magnitude cut (see Sec. l4.2t . no correction is applied. For both fit- 
ting methods, the slopes are also corrected for the (small) effect 
of correlated errors on effective parameters (Sec. I4.3t . using the 
correction factors in Tab. [3] For the magnitude-selected sample of 
ETG, the corrected values of FP slopes are listed in Tabs. [6] and |7] 
for the orthogonal and log ao regression procedures, respectively. 
In Tab.[6]"a" and "b" are the slopes, and "c" and Sr^ are the offset 
and the log dispersion of the FP. Error bars denote la standard 
errors. The quantity s^^ is the intrinsic dispersion of the relation 
along log To. In Tab. |7] "a" and "b" are the slopes, and "c" and 
Sre are the offset and the log dispersion of the FP. Error bars 
denote la standard errors. The quantities "c" and sj,^ are the off- 
set and the scatter of the FP re-measured by fixing the values of 
"a" and "b" in all wavebands. The quantity Sr^ is the amount of 
dispersion in the log rc direction around the plane due to measure- 
ment errors on effective parameters and velocity dispersion, while 
s* ^ is the intrinsic scatter of the FP along the log rc axis. For the 
colour-selected samples, the coefficients are very similar to those 
obtained for the magnitude-selected samples, and are not reported 
here. These tables show how small the statistical uncertainties on 
FP slopes are, amounting to only a few percent in all wavebands. 



Notice that the large number of ETGs makes the NIR FP coeffi- 
cients to have a much better accuracy than any previous study. 

Both the magnitude- and colour-selected samples of ETGs ex- 
hibit very similar trends in Fig. |9] For the log ao fit, we do not see 
any systematic variation of the FP with waveband. From r through 
K, the values of "a" are consistent at less than 2a. In g-band, the 
log ao slope is larger than that in the other bands. The difference 
between g and r-band values of "a" is significant at ~ 3a, after the 
corresponding uncertainties are taken into account|f[ For what con- 
cerns the coefficient "b", all the values are very consistent. On the 
contrary, the orthogonal regression exhibits a clear, though small, 
variation of the slope "a" from g through K. The value of "a" is 
found to vary from ^ 1.38 in g to ~ 1.55 in K, implying a 12% 
variation across the grizYJHK wavebands. The coefficient "b" 
does not change with waveband. We analysed if these results can 
be affected by the (small) contamination of the SPIDER sample 
from early-type spirals. In paper I, we showed that the contami- 
nation from such systems is expected to be ~ 13%. We also de- 
fined a subsample of ETGs with a lower contamination of ~ 5%. 
Fig. [To] compares the FP slopes of the magnitude-selected sample 
with those obtained by selecting only galaxies in the lower contam- 
ination subsample. The values of "a" are fully consistent between 
the two cases in all wavebands, while there is only a marginally 
significant (~2 a) difference in "b" . 

The values of "a" and "b" in Tab. [6] can be compared with 
those obtained from previous studies using the orthogonal fitting 
procedure. The r-band value of "a" is consistent, at 2a, with that 
of a = 1.49 ± 0.05 found by BER03b, and with the value of 
a ~ 1.434 reported bv lHvde & Bemardil ( l2009h . The value of "a" 
is larger, at the 2a level, than that of a = 1.24 ± 0.07 found 



^ To estimate the significance level, we add in quadrature the errors on 
"a" for the two wavebands, assuming they can be treated as independent 
uncertainties. 
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Table 6. Coefficients of the FP in grizY J H K from the orthogonal fit for the magnitude-selected sample of ETGs. 



band a b c Sr^ s* 



9 


1.384 ±0.024 


0.315 ±0.001 


-9.164 ±0.079 


0.125 ±0.002 


0.095 ±0.003 


r 


1.390 ±0.018 


0.314 ±0.001 


-8.867 ±0.058 


0.112 ±0.002 


0.082 ±0.002 


i 


1.426 ±0.016 


0.312 ±0.001 


-8.789 ±0.053 


0.110 ±0.002 


0.079 ± 0.002 


z 


1.418 ±0.021 


0.317 ±0.001 


-8.771 ±0.072 


0.111 ±0.002 


0.079 ± 0.003 


Y 


1.467 ±0.019 


0.314 ±0.001 


-8.557 ±0.058 


0.107 ±0.002 


0.081 ±0.002 


J 


1.530 ±0.017 


0.318 ±0.001 


-8.600 ± 0.060 


0.111 ±0.001 


0.083 ± 0.002 


H 


1.560 ±0.021 


0.318 ±0.002 


-8.447 ±0.077 


0.117 ±0.002 


0.087 ± 0.003 


K 


1.552 ±0.021 


0.316 ±0.002 


-8.270 ±0.076 


0.118 ±0.002 


0.089 ± 0.002 



Table 7. Coefficients of the FP in grizY J H K from the log (tq fit for the magnitude-selected sample of ETGs. 


band 


a 


b 


c 




c' 








9 


1.615 ±0.032 


0.297 ±0.002 


-9.275 ±0.095 


0.135 ±0.002 


-9.080 ±0.002 


0.128 ±0.002 


0.080 ±0.001 


0.100 ±0.002 


r 


1.476 ±0.029 


0.298 ±0.002 


-8.726 ±0.083 


0.112 ±0.001 


-8.813 ±0.002 


0.115 ±0.001 


0.075 ±0.001 


0.087 ±0.002 


i 


1.456 ±0.027 


0.296 ±0.002 


-8.517 ±0.074 


0.107 ±0.001 


-8.694 ±0.002 


0.111 ±0.001 


0.075 ±0.001 


0.082 ±0.002 


z 


1.445 ±0.026 


0.299 ±0.002 


-8.477 ±0.073 


0.104 ±0.001 


-8.605 ±0.002 


0.108 ±0.001 


0.075 ±0.001 


0.078 ± 0.002 


Y 


1.435 ±0.025 


0.297 ±0.002 


-8.164 ±0.073 


0.099 ±0.001 


-8.353 ±0.002 


0.105 ±0.001 


0.066 ±0.001 


0.081 ±0.002 


J 


1.508 ±0.028 


0.305 ± 0.002 


-8.308 ±0.085 


0.103 ±0.001 


-8.195 ±0.002 


0.102 ±0.001 


0.062 ±0.001 


0.081 ±0.002 


H 


1.474 ±0.025 


0.302 ± 0.002 


-7.966 ± 0.074 


0.105 ±0.001 


-7.991 ±0.002 


0.106 ±0.001 


0.068 ±0.001 


0.082 ±0.002 


K 


1.484 ±0.023 


0.300 ± 0.002 


-7.844 ± 0.072 


0.106 ±0.001 


-7.872 ±0.002 


0.107 ±0.001 


0.067 ±0.001 


0.082 ±0.002 
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Figure 10. Effect of limiting the analysis to the sample of ETGs with lower 
contamination from galaxies with residual disc-like morphological features 
(see paper 1). Filled circles and solid ellipses refer to the results of the 
orthogonal fitting procedure for the magnitude-limited sample of ETGs. 
Dashed ellipses and empty circles are those obtained for the sample with 
lower contamination (see the text). Ellipses denote 2 u error contours. The 
FP coefficients turn out to be consistent between the two cases, from g 
through K. 

by JFK96. As noticed by BER03b, the origin of such difference 
is still not understood, although one may notice that it further re- 
duces when considering the value of "a" = 1.31 ± 0.07 found 
from JFK96 for ETGs in the Coma cluster. For what concerns the 
coefficient "h" of the FP, its value in r band (~ 0.314) is con- 



sistent with that of 0. 328 ± 0.008 fouiid by J FK96, and with the 
value of ~ 0.316 from lHvde & BemardilEoOgh . On the other hand, 
BER03b report a somewhat lower value of "6" = 0.300 ± 0.004. 
For what concerns the NIR FP, t he values of the slopes can be com- 
pared with those obtained from lPahre. de Carvalho. & Djorgovskil 
ll998b.) . who found "a" = 1.53 ± 0.08 and "6" = 0.316 ±0.012, 
still very consistent with our findings. The finding that "h" does 
not change with waveband is in full a greement with what already 
suggested by Pahre, de Car valho. & D iorgovski ( 1998b). 

The above results on the waveband dependence of the FP ex- 
tend the findings of LBM08, who derived the FP in the r and 
K bands for 1, 400 ETGs selected with similar criteria as in the 
present study. For the orthogonal fit, LBM08 obtain "a" — 1.42 ± 
0.05 and "b" = 0.305 ± 0.003 in r band, and "a" = 1.53 ± 0.04 
and "6" = 0.308 ± 0.003 in K band. The values of "a" are fully 
consistent with those reported in Tabs. |6] and [T] while the values 
of "6" are smaller, at 2.5 ct, than those we find here. This (small) 
difference is likely explained by the different correction procedure 
adopted here with respect to that of LBM08. In agreement with 
LBM08, we find that, when considering the log ao fit, one does 
not see any significant variation of FP slopes with waveband. When 
comparing the orthogonal fit results in r and K bands, LBM08 
found a variation of only 8 ± 4% (see the values reported in their 
table 1). Here, considering the r and K band values of "a" in 
Tab.|6l we find a variation of 11±2%. The variation is even smaller, 
amounting to ^ 8.5%, when considering the colour-selected sam- 
ples of ETGs. Both values, are consistent, within the uncertainties, 
with those found by LBM08. 



7.2 Dependence on velocity dispersion estimates and 
magnitude range 

As described in paper I, two alternative velocity dispersion esti- 
mates are available for the entire sample of ETGs, those retrieved 
from SDSS-DR6 and the new values we have measured by means 
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Figure 11. Effect of changing the method to derive velocity dispersions 
on FP slopes. Filled and empty symbols are the values of FP slopes ob- 
tained by using SDSS-DR6 and STARLIGHT (tq values, respectively, for 
the magnitude-selected sample of ETGs. Solid and dashed ellipses plot 2 cr 
confidence contours for the samples with SDSS-DR6 and STARLIGHT 
velocity dispersions, respectively. Different colours mark different wave- 
bands, as in Fig. [9] Notice that scales and labelling are the same as in 
Fig.M 



of STARLIGHT jCid Feinandes et alj|2005h . Fig.[TT]compares the 
FP slopes we derive in the different wavebands when using ei- 
ther one or the other set of ao values. Although we find a good 
agreement among STARLIGHT and SDSS-DR6 ctq values (see pa- 
per I), the FP slopes slightly change when using either one or the 
other source of ao's. In particular, the value of "a" is sistemati- 
cally smaller for STARLIGHT with respect to SDSS-DR6. Aver- 
aging over all the wavebands, the difference amounts to ~ —9%. 
We notice that the r-band value of a — 1.26 ± 0.03 from the 
STARLIGHT ao's matches exactly the value of "a" obtained by 
JFK96 (see Sec. |7]l, implying that the method to measure the ao 
might be one main driver of the differences in FP coefficients be- 
tween BER03b and JFK96. Notice also that the value of "b" is 
essentially independent of the velocity dispersion estimates. In or- 
der to analyze if the waveband dependence of the FP is sensitive 
to the magnitude range where ETGs are selected, we proceed as 
follows. First, we select all the ETGs in the SPIDER sample, with 
photometry available in grizYJHK and reduced smaller than 
3. This selection is the same as for the magnitude-selected sample 
of ETGs (Sec. |2]l, but without applying any magnitude cut in r- 
band. The sample consists of 4, 981 galaxies. Fig.[T2]compares the 
FP slopes of the magnitude-selected sample of ETGs with those 
obtained for the entire sample. 

The slopes of the FP are fully consistent in all wavebands be- 
tween the two cases. We also define two subsamples consisting of 
all the ETGs with available photometry in grizY J H K and r-band 
magnitude brighter that Mr= -21 and °-'^Hlr= -21.5, re- 
spectively. We exclude galaxies whose Sersic model fit gives an 



high X value (> 3). These 



'Mr 



-21 and "■"'M^= -21.5 
subsamples include N = 411 and N — 2, 091 galaxies, respec- 
tively. Fig.llSlcompares the slopes of the FP obtained for these two 
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Figure 12. Same as Fig.[To]but showing the effect of using the entire SPI- 
DER sample, rather than the magnitude-selected sample of ETGs, on the 
waveband dependence of the FP. Notice that scales and labelling are the 
same as in Fig. 1101 
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Figure 13. Effect of changing the magnitude limit on the waveband depen- 
dence of the FP. The upper and lower panels plot as open circles, connected 
by dashed segments, the FP slopes "a" and "h" obtained when apply- 
ing different magnitude limits of ^'^"^ Mr=— 21 and —21.5, respectively. 
In both panels, the results, obtained for the entire magnitude-selected sam- 
ple, are also shown, for comparison, as filled circles, connected by solid 
segments. Different colours mark different wavebands, as in Fig.|9] 
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samples, by the orthogonal fitting procedure, with those obtained 
for ETGs in the magnitude range of °-°'^Mrfk -20.55. In order to 
allow a direct comparison of the amounts of variation with wave- 
band, for both subsamples, the best-fitted values of "a" and "h" 
are rescaled to match the values of "a" and "h" in the r-band for 
the magnitude-selected sample. The figure clearly shows that the 
waveband dependence of the FP is essentially the same regardless 
of the magnitude range. For ""^A/,.^ —21, the variation of "a" 
is smaller, but consistent within the errors, with that obtained for 
°-^Hlr4, -20.55. For Mr^ -21.5, the trend of "a" vs. "h" 
matches very well that obtained for the entire sample. In all cases, 
the log (JO coefficient increases from g through K, while the value 
of "b" is independent of waveband. 

7.3 Dependence on galaxy parameters 

The FP relation and its dependence on waveband might change 
when selecting samples of ETGs with different properties. To an- 
alyze this aspect, we split the magnitude- selected sample accord- 
ing to the value of different galaxy parameters, i.e. the axis ratio, 
b/a, the Sersic index, n, the r — K colour index, and the average 
discy/boxiness parameter, 04. We utilize the values of b/a, and n, 
in the r band, while for 04, we adopt its median value among the 
gri wavebands (see paper I). The r ~ K colour is computed from 
2DPH0T total magnitudes. 

Fig.Elplots the slope "a" of the FP as a function of "b". The 
slope's values are those obtained from the orthogonal fit, applying 
the same correction factors as for the entire sample (Sec. 17. It . Each 
panel corresponds to a given parameter p: a^, b/a, n, and r--K. For 
each parameter, the magnitude-selected sample of ETGs is splitted 
in two subsamples, having values of p either lower or higher than a 
given cut value, pc- For p = b/a, n, and r ~ K,we set pc equal to 
the median value of the distribution of p values. The median values 
arepc ~ 0.699, 6.0 and 3.0 for b/a, n, and r — A', respectively. For 
a4, we divide the sample into discy (04 > 0) and boxy (04 < 0) 
galaxies. Notice that, for a given parameter, galaxies in the two sub- 
samples can populate different regions of the log rc-< n >c-log uo 
space. For instance, because of the luminosity-size relation and the 
KR, galaxies with higher Sersic index are brighter and tend to have 
higher values of < /i >e. This geometric difference might produce 
spurious differences in FP coefficients. A trivial example of this 
geometric effect is the magnitude selection: the bias on FP coeffi- 
cients changes for samples of ETGs spanning different luminosity 
ranges (Sec. l4.2l (. In order to minimize any geometric difference, 
for a given parameter p, we extract the two subsamples of ETGs by 
constraining their distributions in magnitude and < /i >c to be the 
same (see App.|B]for details). 

Fig.ll4lshows that the waveband dependence of the FP is sim- 
ilar for all subsamples, i.e. the value of "b" tends to be constant 
while the coefficient "a" increases by ~ 15% from g through K. 
However, the FP slopes change significantly for samples of ETGs 
with different properties. The differences can be summarized as 
follows. 

- Galaxies with higher n have a lower value of "b"; the value of 
"a" in the MR is smaller for the subsample with n > 6, while in 
the optical both subsamples have consistent "a" . 

- The FP of round galaxies (higher b/a) is more tilted (smaller 
"a") than that of galaxies with low b/a. The difference is more 
pronounced in the NIR than in the optical. 

- For a4 and r — K, one can notice a different behaviour. In the 
NIR, the FP slopes of the two subsamples are fully consistent, while 



in the optical, there is a detectable difference in the coefficient "b" . 
Boxy and blue (i.e. r ~ K <?,) galaxies tend to have lower "b". 

We remark that all these trends remain essentially unchanged when 
replacing SDSS-DR6 with STARLIGHT velocity dispersions, with 
the exception that "a" is slightly lower for STARLIGHT relative 
to SDSS (70 values (see Sec. l7.2t . 

Fig.[T5]shows the FP slopes obtained for different subsamples 
as in Fig. [141 but without imposing the constraint that, for a given 
quantity, the two subsamples consist of galaxies with the same dis- 
tributions in magnitude and < /i >c. No difference would have 
been detected with respect to n and b/a, while a (spurious) differ- 
ence in the NIR value of "a" between red and blue galaxies would 
have been found. The comparison of Fig. [15] and Fig. [14] proves 
that accounting for purely geometric differences in the space of 
FP parameters is of paramount importance to correctly analyze the 
scaling relations of different galaxy samples. 



8 THE EDGE- AND FACE-ON PROJECTIONS OF THE 
FP 

So far, we have analysed the waveband dependence of the FP, and 
that of the FJ and KR. Since the FJ and KR are projections of the 
FP, we expect their waveband dependence to be connected to that of 
the distribution of galaxies in the FP. We establish this connection 
by analysing the edge- and face-on projections of the FP. 

Fig.[T6]presents the so-called short edge-on projection of the 
FP, from g through K, namely the combination of effective struc- 
tural parameters, log re — 6 < /i >e, as a function of logao. 
This corresponds to the FP along the shortest axis, whose slope 
is equal to its logtro coefficient, "a". Each panel in Fig. [16] for a 
given passband, shows the FP obtained from the orthogonal fitting 
method (solid line), as well as the r-band fitted FP (dashed line). 
From these plots we can see the increasing of "a" from the optical 
through the NIR. Comparing the solid and dashed lines we see that 
the increasing is quite small (see Sec. l7.lt . The observed scatter in 
the edge-on projection decreases from the optical through the NIR, 
as it can be attested from the values of the FP log dispersion, 
Sre > reported in Tab. [70. 

In order to represent t he FP face- on projection, we follow the 
same formalism as in Guz man. Lucev . & Bower ( 1993) (hereafter 
GLB93). We project the FP into a plane defined by two orthogonal 
directions, one of which is perpendicular to the logrc axis. The 
axes of the projection are: 

X' = (xologre + 6' log < J>e +a logao) /■\/xo x (1 + ig) 
Y' = (alog<J>e -&'logao)/V^, (7) 

where b' = —6 x 2.5, xo = + {b'Y , and < J >5 is the mean sur- 
face brightness in flux units, with < /i >e= —2.5 log <I>e- From 
the FP equation and Eq.[6] it follows that X' is simply proportional 
to logro. Fig. [TT] shows the distribution of ETGs on the face-on 
projection of the FP in r-band, together with the log rc, < p >o, 
and log (To, directions, as well as the directions of increasing total 
magnitude, A'l AG, and logarithmic luminosity, log L, on the face- 
on FP. The dashed lines in the plot illustrate the oq and magnitude 

We actually refer to the .scatter values reported in Tab.|2](and not those 
reported in Tab. |6) as for the log ctq fitting procedure the FP slopes do not 
change significantly with waveband, allowing a meaningful comparison of 
the Sr^ values from g through K. 
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Figure 14. FP slopes, from g thi'ough K, for different subsamples of ETGs. Each panel shows the FP slopes obtained by splitting the magnitude-selected 
sample of ETGs in two subsamples, according to the Sersic index (panel a), the axis ratio (panel b), the isophotal parameter, 04 (panel c), and the r — K 
colour index (panel d). For each quantity, the two bins are defined as shown in the lower-right comer of the coiresponding panel. The slope's values are plotted 
with different symbols and are connected through different line types for the two bins, as shown in lower-right comer of each plot. Different wavebands are 
represented with different colours, as in Fig. [9] Ellipses denote Ict confidence contours for "a" and "6". To make the plot more clear, only the ellipses in g- 
and /('—bands are shown. Notice that scales and labelling of each panel are the same as in Fig. II 01 



selection limits of the sample (Sec.[2ll. As already noticed in previ- 
ous studies (e.g. lBender. Burstein. & Fabeill992f : GLB93; JFK96), 
ETGs are confined in a small region of the face-on projection, only 
partly due to selection effects. Galaxies populate a diamond-shaped 
region, limited at low X' by the magnitude limit of the sample, 
Mr,iim = —20.55, and at high X' , by the bright-end knee of the 
galaxy luminosity function, i.e. the fact that there are no galaxies 
brighter than a magnitude threshold of about Mr,um — 4. Notice 
that the ctq selections (see paper I and Sec. ^ do not affect the 
shape of the distribution in the face-on projection, as all galaxies 
lie well within the region defined by these additional cuts (dashed 
green lines in the figure). 

Since the logrc and < /i >c directions form almost a 90° 
angle on the FP (see the blue and magenta arrows in the upper-right 
of Fig.llVt. the KR is essentially reflecting the face-on distribution, 
as already noticed by GLB93. In order to establish this connection, 
we perform an orthogonal least-squares fit of the diamond-shaped 
region, accounting for the magnitude selection in the X'-Y' plane 
by the MLSO fitting procedure (see Sec. [5). The relation is: 



Y' = const + A' X X' , 



(8) 



obtain a best-fitting value of A' = -1.08 ± 0.01. Since the log rc 
and < >c directions are approximately orthogonal, the fitted line 
is very similar to what we would obtain by binning the data with 
respect to log and take the median values of X' and Y' in each 
of those bins. The result of this binning procedure is shown by the 
magenta circles in Fig. [17] The magenta line is the best fit of the 
binned data-points, with a slope of -1.01±0.02, very close to the 
MLSO fit reported above. The 2 a scatter of the MLSO fit, along 
the X' , is displayed by a segment in the lower-left of Fig.[T7l with 
the shorter segment corresponding to the 2 a log rc dispersion of 
the FP seen edge-on (Tab.[6]l. The scatter around the edge-on FP is 
about twice smaller than that of the face-on FP as already noticed 
by GLB93 implying that the FP is more like a band rather than 
a plane, in the logro, < p >c, logcro space. We can use Eq. [8]to 
eliminate log ctq from the FP equation (Eq.[T). This leads to a linear 
relation between < fi>c and logro, similar to Eq. |2] i.e. the KR, 
whose expected slope is: 



P2 



ax A X . 

Xo QAxo 



(9) 



where A' is the slope, and const is an offset. For the r band, we 



Inserting the r-band value of the FP slope from Tab.[6]and the best- 
fitting value of A' in this equation, we obtain p'2 = 3.56 ± 0.03, in 
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Figure 15. The same as in Fig. 1141 but without matching the distributions in magnitude and mean surface brightness of subsamples in the two bins of a given 
quantity. 
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Figure 16. Short edge-on projection of the FP (see Sec. [s), where the photometric quantity entering the FP, logre — b < fi >e, is plotted against the 
spectroscopic quantity logtro. Different panels correspond to different passbands, from g (upper-left) through K (lower-right), as indicated in the upper-left 
comer of each plot. In the short edge-on projection, the FP projects into a line, having a slope equal to the FP coefficient "a". For each panel, this projection 
is shown by the solid light-grey line. In order to emphasize the waveband dependence of "a", in each panel we plot as a reference, with a dashed dark-grey 
line, the r-band FP projection. Notice that the value of "6", defining the y-axis variable changes among different panels, according to the values reported in 
Tab.[6] Notice that for a more direct comparison of the FP projection in different wavebands, the lengths of the x- and y-axes are the same for all panels. 
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Figure 17. Face-on projection of the FP in r band. Tlie projection is such that the x-axis variable, X', is proportional to logrc, while the y-axis variable, Y', 
is proportional to — < fi >e and logco- The aiTows in the upper-right corner of the plot denote the directions where the quantities, MAG, logcro, logL, 
<fj,>c, and log r'c increase, where AIAG is total galaxy magnitude, and log L is logarithmic luminosity in r band. The size of the arrows amounts to 0.5 dex, 
1 mag arcsec~^, 3 mag, and 0.5 dex for log (tq, < >c, MAG, and log r'c, respectively. The red dashed lines correspond to the r-band magnitude hmit, 
Mr,iim = —20.55, and a bright-end limit four magnitudes brighter than M,, ii^, as shown by the corresponding labels. The green dashed lines correspond 
to the log (TO lower and upper selection limits of 70 and 420A:ms~^. The solid cyan line is the MLSO best-fitting relation to the data (see the text). The 
thick solid black line marks the exclusion zone, Y' = — 0.59X' + const, as originally defined bv lBender . Burstein, & Fabei 1992. The intercept of this line 
has been arbitrarily normalized to mark the upper envelope of the face-on projection of the FP. The presence of a few points above the line is likely due to 
measurement errors on FP variables. The magenta and red circles are obtained by binning the data with respect to log and log L, respectively, computing 
the median values of X' and Y' in each bin. The magenta and red solid lines are the best-fitting fines to the binned data-points. The size of the 2 a scatter 
around the fit of the face-on projection (cyan line) is given by the long segment in the lower-left comer of the plot. Notice how this segment is about twice 
larger than the corresponding log scatter of the FP, given by the short segment at lower-left. 



good agreement with the KR r-band slope, p2 ~ 3.55, obtained by 
theMLSB fit (see Tab. EJ. 

As far as the FJ relation, we notice that the log L and log ctq 
axes form a small angle on the FP, and are almost orthogonal to 
the direction of the long diagonal of the diamond-shaped region. 
As a consequence, the best-fitting line of the face-on distribution 
cannot be directly connected to the FJ, as we do for the KR. In 
fact, the FJ relation almost coincides with the (short) edge-on pro- 
jection of the FP (e.q. GLB93). In order to overcome this problem 
and relate the face-on distribution and the FJ relation, we bin the 
data with respect to log L and then compute the median values of 
X' and Y' in each bin. This binning procedure allows us to look 



at the distribution of galaxies in different luminosity bins, in the 
same way as for the FJ relation. The log L-binned points are plot- 
ted as red circles in Fig. [TT] The corresponding linear best-fitting 
is displayed as a red line. The slope value of the red line amounts 
to A'l^ = -0.72 ± 0.03. Notice how the MLSO fit and the red line 
differ significantly. Replacing A' with A'j^ in Eq.[8] and combining 
the resulting equation with the FP, we obtain a linear relation be- 
tween log L and log ctq, similar to the FJ equation (Eq.[5j, with an 
expected slope of 

y ^ a - b' X A'l X ^1 + XQ ^^^^ 
^ a X X \/T+~Eo + b' + 2x0 
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Table 8. Slopes of the face-on projection of the FP, A' and A'j^, and predicted 
slopes of the FJ and KR, A ^ and p'^ . 



band 


A' if it) 


^1 


V'2 


K 


9 


-1.05 ±0.01 


-0.77 ±0.03 


3.48 ±0.02 


0.11 ±0.01 


r 


-1.08 ±0.01 


-0.72 ±0.03 


3.56 ±0.03 


0.14 ±0.01 


i 


-1.07 ±0.01 


-0.69 ±0.03 


3.51 ±0.03 


0.15 ±0.01 


z 


-1.09 ±0.01 


-0.78 ±0.04 


3.54 ± 0.03 


0.10 ±0.02 


Y 


-1.14 ±0.02 


-0.67 ±0.04 


3.64 ± 0.05 


0.16 ±0.02 


J 


-1.18 ±0.02 


-0.68 ±0.03 


3.66 ± 0.06 


0.14 ±0.02 


H 


-1.17 ±0.02 


-0.73 ±0.04 


3.67 ±0.06 


0.13 ±0.01 


K 


-1.24 ±0.01 


-0.70 ±0.04 


3.83 ± 0.04 


0.14 ±0.01 



Inserting the value of A'j^ and the FP coefficients in this equation 
we obtain A'l = 0.14 ± 0.01, very close to the measured slope 
of the FJ relation in r band (Ai = 0.19 ± 0.02, see Tab.©. The 
difference between Ai and A'l does not reflect any inconsistency 
in the data, but just the fact that the log L and log ctq directions 
form a small angle on the FP, and hence it is not straightforward 
to connect the distribution of galaxies on the face-on projection to 
that on the log L-log (tq plane. Eq.[T0]is used here as an empirical 
tool to analyse the dependence of the FJ relation on waveband (see 
below). 

Fig.[T8]shows the face-on projections of the FP from g through 
K. For each band, we have performed an MLSO fit of the data, as 
well as a log L-binned fit, in the same way as we do in Fig. [T7] 
For each band, the corresponding slopes. A' and A^, are reported 
in Tab. [8] together with the predicted slopes of the KR, p'2, and FJ 
relation, A'l, from Eqs.|9land|10| respectively. From Tab.[8]we see 
that the slope of the KR is expected to increase with waveband, in 
agreement with what we measure (Sec.O. This can be seen directly 
from Eq.|9l as "b" does not change significantly with waveband, 
and the same holds for the term ^/l-^- Xq / Xq . It follows that the 
waveband dependence of the KR slope is driven by the term —a x 
A' (second term of Eq. |9ll. From the values of A' in Tab. [S] we 
see that —A' increases with waveband, i.e. the MLSO fitted line 
steepens with waveband, in the same way as "a" does (Sec. 17. il l. 
Therefore, the variation of the KR from g through K is connected 
to the variation of both the slope, "a", and the face-on projection 
of the FP with waveband. The steepening of the MLSO fit from g 
through K can be explained by the variation of optical to NIR radii 
along the ETG's sequence (Sec. |5]l, and considering the fact that 
X' is essentially proportional to log r^. 

For what concerns the FJ relation, the slope listed in Tab. [8] 
does not change with wavelength, which is consistent with the re- 
sults presented in Sec.|6] Eq. [TO] explains the reason for this be- 
haviour. First, we notice that the FP slope, "a", appears both in the 
upper and lower part of the second term of Eq. [TO] Therefore, the 
waveband dependence of "a" does not affect the A'^ value. More- 
over, from Tab. [H we see that the log L-binned slope of the FP 
face-on distribution is independent of waveband, making the value 
of A'l constant from g through K. In other words, the face-on distri- 
bution of the FP changes with waveband in a complex way, so that 
the long diagonal of the diamond-shape region steepens with wave- 
band, while the log L-binned envelope of the distribution does not 
change with waveband. The former effect, together with the wave- 
band variation of the FP coefficient, "a", determines a dependence 
of the KR on waveband, while the latter is consistent with the FJ 
relation not changing from g through K. 



9 STELLAR POPULATIONS ALONG THE FP 

Under the homology assumption, one can combine the FP relation 
(Eq.[TJ with the virial theorem 



2 M 
aooc- 



</>e r-e 



(11) 



and parametrize the mass-to-light ratio, as a function of two 
variables out of M, L, ctq, r^, and < I (Djorgovski, de Car- 
valho, & Han 1988). Here, we denote as < J >e the mean surface 
brightness within re in flux units. In order to analyze how stellar 
population parameters vary along the sequence of ETGs, it is con- 
venient to parametrize such sequence by means of variables that 
are independent of stellar population parameters. To this effect, we 
consider the quantities M and cro, and write 



— OC 

Li 



(12) 



where the index x runs over all the available wavebands (x = 
grizYJHK). Using Eq. [T]and Eq. [TT] one obtains the following 
expressions for and fix'- 



= 4-0.4 



(^) 



(13) 
(14) 



04 

where ax and bx are the values of the logao and < ^i>c slopes 
of the FP in the waveband x. These equations imply that, at fixed 
(Jo the variation of the M / L with mass is completely determined 
by the coefficient "b" of the FP. On the contrary, at fixed AI, the 
variation of M/L with velocity dispersion is determined by both 
the values of "a" and "b" . Hence, the result that "b" does not 
change from g through K (see Sec. |7j implies that, at fixed ao, 
the change of M/L with mass is independent of waveband. On the 
contrary, the dependence of A'l/L with ctq (at fixed AI) changes 
from g through K. This is shown in Fig. [19] where we plot the val- 
ues of a and /3 in the grizYJHK wavebands. For each band, we 
calculate a and /3 from Eqs.[T3]and[T4l using the FP coefficients ob- 
tained by the orthogonal fitting procedure (see Tab.|6]l. The reason 
for adopting the FP slope's values from the orthogonal regression is 
discussed in Sec. 9.1. As expected, the value of /3 is constant, while 
a increases from g through K. Although the variation of "a" from 
g through K amounts to only ~ 12%, the corresponding increase 
in the a value is significant, amounting to ~ 70%. 

One can also notice that the values of a and /3 have opposite 
sign. Since the value of a is negative, at a given mass, the AI/L is 
a decreasing function of ao. On the contrary, for fixed ao, the AI/L 
increases with AI. In order to characterise the overall variation of 
AI/L along the ETG's sequence, parametrized in terms of galaxy 
mass, we have to project Eq.[T2] i.e. the FP itself, into the M/L-M 
plane. To this effect, we can take advantage of a specific projection 
of the FP, such as the FJ relation, i.e. the fact that luminosity is 
proportional to ao. In this approach, the FJ relation is not providing 
any extra information wrt the FP itself, but is used as an empirical 
tool to project Eq.[T2]into an AI/L vs. AI power-law[fl Using Eq.|5] 
to replace ao with L in Eq.[T2] for a given waveband X, we obtain: 



^ A different approach would be that of measuring directly dynamical 
mass from the data, by means of the virial theorem. This approach (see 
e.g. JFK96) relies on a given galaxy model to translate cro and into M, 
and hence implies several assumptions about, for instance, the dark-matter 
component of ETGs. This analysis is currently under way for the SPIDER 
sample, and will be presented in a forthcoming contribution. For the present 



SPIDER - II. The FP in grizYJHK 1 9 




Figure 18. The same as Fig.[T7]for the grizY J H K wavebands, from left to right, and top to bottom. For better displaying the plots, the internal labels of 
Fig. ll7l are not shown. Panel (b) is the same as Fig. ll7l and is repeated to allow a direct comparison with the other panels (wavebands). Notice that for a more 
direct comparison of the face-on projection in different wavebands, the lengths of the x- and y-axes are the same in all panels. 



M/L oc 



(15) 



M/L oc M 



(18) 



where 



7^ 



(1 



A) 



(16) 



and A = 0.198 ± 0.007 is the average slope of the FJ relation 
(Sec.[6ll. We point out that using the values of A we have measured 
for each waveband (see Tab. 5) rather than the average A value 
would not change at all the results presented here. The values of 7, 
derived with the above procedure, are reported in Tab.|9](col. 2) for 
each waveband. The 7 has a positive value in all wavebands, and 
tends to slightly decrease, by 20 ± 5%, from g to K. This variation 
can be interpreted as a change of stellar population properties along 
the sequence of ETGs. To this effect, following the same approach 
of LBM08, we assume that also the stellar mass-to-light ratio of 
ETGs, M*/L, is a power-law of M: 



Eq.[T2]can then be written as 



(17) 



study, we adopt a model-independent approach, using only the information 
provided by the FP. 



where 7' = 7a; — 7^ defines how the ratio of stellar to total mass 
changes along the mass sequence of ETGs, M,/M oc and 
thus it is assumed to be independent of waveband. Introducing the 
parameter / = '^k/ik, which defines the fraction of the K-band 
slope of the M/L vs. M relation due to stellar population effects, 
we obtain the following system of equations: 



1-/ 



IK +lx ^ "fx- 



(19) 



We note that / can vary between and 1. For f — 0, the K-band 
tilt is independent of stellar populations (7J- — 0), while for f — 1 
the tilt is entirely explained by stellar population effects (7]^ — 7x 
and 7' — 0). The quantities 7* depend on how stellar population 
properties change along the mass sequence of ETGs. Considering 
only the age, t, and the metallicity, Z, one can write: 



. _ Sjlog M,/L) _ Sjlogt) 
" (5(log M) ^ ' 5(log M) 



where the quantities S{logt) and 5{\ogZ) are the logarithmic 
differences of age and metallicity between more and less mas- 
sive galaxies (per decade in mass), while Ct^ = 



81ogAi,/£^ 

aiogt 



and 
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Figure 19. Slopes, a and fi, of the power-law of -j- as a function of ctq 
and M in the grizY J H K wavebands. Solid and dotted ellipses denote Icr 
and 2(r confidence contours, respectively, as implied by the uncertainties on 
FP coefficients (see Tab.[6). 



Cz^ = ^ ^°eiog,z'^'° '"^^ partial logarithmic derivatives of M^/L 
(in the waveband x) with respect to t and Z. Deriving the coeffi- 
cients 'yx from the slope's values of the FP in the different wave- 
bands (Eqs. [13] [T4I and I16t . and inserting the expression of 7^ 
from Eq. [20] into Eq. [191 we obtain a system of eight equations, 
one for each of the grizY J H K wavebands, in the three unknowns 
(5(log t), 5(log Z), and /. We solved this system by minimizing the 
sum of relative residuals: 



X = 



E 



(1 - /)/./ -iK+ll-l^ 



(21) 



We estimated the quantities Ct.x and cz,x using simple stellar pop- 
ulation (SSP) models from diffe rent sources: iBruzual & Chariot 
( I2OO3I) (BC03). lMarastonl j2005l) (M05), and Chariot and Bruzual 
(2010, in preparation; CBIO). These models are based on different 
synthesis techniques and have different IMFs. The M05 model uses 
the fuel consumption approach instead of the isochronal synthesis 
of BC03 a nd CB 10. The CB 10 cod e implements a new AGB phase 
treatment jMarigo & Girardill2007l) . The IMFs are: Scalo (BC03), 
Chabrier (M05), and Salpeter (CBIO). Moreover, we also used a 
composite stellar population model from BC03 having exponential 
star formation rate (SFR) with e-folding time of r = 1 Gyr (here- 
after _BC03t=i). The models are folded with the grizYJHK 
throughput curves, and the M / L values computed for different val- 
ues of t and Z. In order to evaluate the impact of changing t and Z, 
we considered three cases, with (i) an age of 9 Gyr and solar met- 
alicity, (ii) an older age of 12 Gyr and solar metallicity, and (iii) an 
age of 12 Gyr and super-solar metallicity (Z = I.^Zq). The min- 
imization was performed 2000 times for each kind of model, and 
for each combination of t and Z values, shifting each time the FP 
coefficients according to the corresponding (correlated) uncertain- 
ties. For each iteration, we found that all the eight equations were 
solved with an accuracy better than 10%. 
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Figure 20. Best-fitting values of 5(log t), 5(log Z) and / values. The left 
panels plot (5(log t) vs. 5(log Z)forthe BC03, CBIO, M05, and BC03t=i 
SSP models (from top to bottom). Different colours denote the different 
combinations of t and Z as shown in the upper-light comer of the top panel. 
For a given panel and colour, the different points correspond to the solutions 
obtained in the 2000 minimisation iterations. The solutions corresponding 
to /< 0.1 are plotted with larger symbols. The light panels plot the corre- 
sponding distributions of / values. To make the plot more clear, small shifts 
(of ±0.05) have been applied to the histograms with different colours. 



Fig.[20|plots 5(logt) vs. (5(log Z), as well as the distribution 
of / values obtained in each case for all the 2000 iterations. For al- 
most all solutions, the / is very close to zero, implying that the tilt 
of the NIR FP is not due to a variation of stellar population proper- 
ties of ETGs with mass. For instance, in the case t = 9Gyr and so- 
lar metallicity, the percentages of solutions with / < 0.05 amounts 
to 94%, 82%, 95%, and 94% for the BC03, CBIO, M05, amd 
BC03t=i models, respectively. Considering only the solutions with 
/ < 0.05, we estimated the mean value of (5(log t) and (5(log Z). 
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Table 9. Slopes of the M/L vs. mass relation, obtained by projecting 
the FP through the FJ relation (col. 2), by using STARLIGHT (rather 
than SDSS) tro 's (col. 3), and fitting the 7 to the FP coefficients rather 
than using the FJ relation (col. 4). 



waveband 




7 






SDSS (JO 's 


STARLIGHT cto's 


a-/3fit 


(1) 


(2) 


(3) 


(4) 


g 


0.224 ± 0.008 


0.251 ±0.009 


0.249 ± 0.008 


r 


0.225 ± 0.006 


0.253 ±0.007 


0.248 ± 0.007 


i 


0.221 ±0.006 


0.247 ±0.006 


0.254 ±0.007 


z 


0.213 ± 0.008 


0.236 ±0.006 


0.233 ± 0.007 


Y 


0.208 ± 0.007 


0.230 ± 0.007 


0.227 ± 0.007 


J 


0.186 ±0.007 


0.202 ±0.008 


0.215 ± 0.008 


H 


0.180 ±0.009 


0.221 ±0.007 


0.208 ± 0.007 


K 


0.186 ±0.009 


0.218 ±0.007 


0.214 ±0.008 



The mean values, and the corresponding uncertainties, are reported 
in Tab. [TO] The uncertainties were estimated by the standard de- 
viation of the (5(logt) and S{\ogZ) values obtained for a given 
model and for a given combination of t and Z. The mean values 
do not depend significantly on either the model or the adopted val- 
ues of t and Z. On average, S{logt) is very close to zero, while 
the (5(log Z) mean value amounts to ~ 0.1 dex. This implies that 
ETGs have syncronous luminosity-weighted ages, with an age vari- 
ation smaller than a few percent per decade in mass, while the 
metallicity variation per mass decade amounts to ~ 23%. These 
results remain unchanged when using STARLIGHT rather than 
SDSS velocity dispersions. Inserting the values of FP coefficients 
and FJ slope as estimated by STARLIGHT (rather than SDSS) ao 's 
in Eq. [I6l the values of the 7's increase on average by ~ 12%, as 
shown by comparing the values in column 3 (STARLIGHT) and 
column 2 (SDSS) of Tab. |9] For instance, the value of 7g changes 
from 0.224 (SDSS) to 0.251 (STARLIGHT), while in K band the 7 
varies from 0.186 (SDSS) to 0.218 (STARLIGHT). Applying the 
procedure described above to estimate /, 5(logt) and 5{logZ), 
we find that, for BC03 SSP models with an age of 9 Gyr and so- 
lar metallicity, the corresponding distribution of / values is still 
strongly peaked around zero, while the mean values of S{logt) and 
5(log Z) amount to ~ 0.008 and ~ 0.100, respectively, fully con- 
sistent with what obtained from SDSS ctq's (<5(logt) ~ 0.013 and 
5(log Z) ~ 0.105, see Tab. I lot. As a further test, we estimated the 
7's by not using the FJ relation. Combining Eq.[T7]with the virial 
theorem, we obtain the equation: 

l-lx 0.4 
log r-c oc 2 X log ao + — < ^ >o, (22) 

1 + 7^ 1 + 7^ 

which reduces to the FP for a^; — 2(1 — 7^)/(l + 7a^) and 
bx ~ 0.4/(1 + 73:). For each waveband, we estimate the by 
minimising the expression: 

= (a.. - 2\-^) /{Sa,f+ (b^ - /{5b^)\ (23) 

where and bx are the FP coefficients from Tab. |6l and Sa^ and 
5bx are the errors on a^; and b^. The corresponding values of 7^; 
are reported in Tab.|9](col. 4). On average, the values of 7^ tend to 
increase, wrt to those from Eq.[T6l by ~ 13%. Even in this case, 
this variation does not impact at all the above conclusions, i.e. the 
/ is zero, while the mean values of 5(\ogt) and 5{\og Z) amount 
to about 0.01 and 0.1 dex, respectively. 



Table 10. Age and metallicity differences per decade of galaxy mass. 



MODEL t{Gyr) Z/Zq (5(logt) 5(\ogZ) 





12 


1 


0.013 ±0.021 


0, 


104 


± 


0, 


,026 


Bcm 


9 


1 


0.013 ±0.017 


0, 


,105 


± 


0, 


,025 


Bcm 


12 


1.5 


0.004 ± 0.001 


0. 


,106 


± 


0, 


,019 


CBW 


12 


1 


0.003 ±0.021 


0. 


,121 


± 


0, 


,022 


CBW 


9 


1 


0.005 ±0.025 


0, 


121 


± 


0, 


,026 


CBW 


12 


1.5 


-0.019 ± 0.004 


0, 


150 


± 


0, 


,027 


A/05 


12 


1 


0.008 ±0.018 


0, 


112 


± 


0, 


,025 


M05 


9 


1 


0.012 ±0.023 


0, 


108 


± 


0, 


,022 


M05 


12 


1.5 


0.005 ±0.001 


0, 


,094 


± 


0, 


,017 


BC03r = l 


12 


1 


0.011 ±0.012 


0, 


,107 


± 


0, 


,023 


BC03r=i 


9 


1 


0.012 ±0.014 


0, 


107 


± 


0, 


,025 


BCOZr = l 


12 


1.5 


0.006 ±0.001 


0, 


100 


± 


0, 


,018 



10 DISCUSSION 

10.1 The fit of the FP in different wavebands 

One of the crucial aspects of the present study is the fitting pro- 
cedure used to obtain the coefficients of the FP and how they are 
affected by different systematic effects. Different fitting techniques 
produce different estimates of FP coefficients, and may lead to er- 
roneous results when co mparing the FP relations obtained with dif - 
ferent samples (LBC00: ISagUa et al.ll200ll : iBernardi et alj|2003bh . 
To avoid this problem, we adopt the same fitting method for all 
the ETG subsamples we analyze. Selection effects and correlated 
errors on effective parameters can be taken into account analyt- 
ically under the assumption t hat the FP variables are normally 
distributed ( iSaglia et alJboOlh . Although IBernardi et all ( I2003H) 
showed that the joint distribution of log Tc, log ctq, and galaxy mag- 
nitude is relatively well described by a multivariate Gaussian, this 
might not necessarely be true when effective parameters are de- 
rived by the Sersic (2DPH0T) rather than de Vaucouleurs (Photo) 
model. These two pipelines yield significant differences in log 
and magnitudes (see paper I). These differences depend on galaxy 
magnitude, and may be partly due to the sky overestimation prob- 
lem affecting the SDSS Photo parameters. We have adopted a non- 
parametric approach, first fitting the FP relation and then correcting 
the slopes for different systematic effects using extensive Monte- 
Carlo simulations. We find that the main source of bias on the FP 
slopes is the magnitude cut. In agreement with iHvde & Bemardil 
('2009), we show that for the orthogonal fit this cut leads to un- 
derestimating the FP coefficients, with the effect becoming negli- 
gible only at faint magnitude limits {Mr ~ —18.5, see Fig. |2]l- 
The effect is negligible when we use the logao fitting method. 
As shown by LBCOO, minimizing the logao residuals leads to a 
log (Jo slope of the FP systematically higher than that obtained by 
other fitting techniques (see also JFK96). We also find that the coef- 
ficient "a" of the FP in the optical (SDSS) wavebands are system- 
atically larger when we use the log ao method compared to results 
obtained with the orthogonal fitting procedure. In r-band the differ- 
ence amounts to ~ 6%. On the other hand, the coefficient "b" turns 
out to be systematically lower, by ~ 5%, for the log ctq method 
(see Tabs. [7] and [6j. Another important result we find is that the 
difference produced by different fitting method depends on wave- 
band (see also LBM08). The FP coefficients do not change with 
the waveband when using the log ao method, while they smoothly 
vary, by ~ 12%, from g through K when using the orthogonal 
method. This can be explained by the fact that the log ctq regres- 
sion minimizes the rms of residuals in the perpendicular direction 
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to the log 7-e-< >c plane, and hence it is less sensitive to differ- 
ences in the distribution of galaxies in that plane, like those among 
effective parameters measured in different wavebands. The prob- 
lem of deriving the best fitting coefficients of correlations among 
astrophysical quantities has been addressed bv llsobe et alj ( fl990h . 
They concluded that, in case one aims to study the underlying func- 
tional relation among the variables, regression procedures treating 
all the variables symmetrically, like the orthogonal method, should 
be adopted. For this reason, we have analyzed the implications of 
the waveband dependence of the FP adopting the results of the or- 
thogonal regression. 

10.2 Variation of roPT /tnir witli galaxy radius 

In the present study, we find that the slope of the KR exhibits a 
small systematic variation with waveband, steepening by ~ 10% 
from g through K. This variation may be explained as the ratio of 
optical to NIR effective radii decreasing for galaxies with larger 
re, namely, while smaller size ETGs have, on average, optical radii 
larger than the NIR ones, the most massive galaxies have r^pj. ~ 
r-j^jj^.ln the assumption that r^p^ /'''nir ^ proxy for the internal 
colour gradient of an ETG, this finding implies that the stellar popu- 
lations of the most massive ETGs have a more homogeneous spatial 
distribution insid e the galaxies, i . e. fl atter radial gradients, than less 
massive systems. ISpolaor et alj found that the relation be- 

tween the internal metallicity gradient and mass in early-type sys- 
tems is bimodal, with a sharp transition at Mb ~ — 19. This magni- 
tude corresponds approximately to the lower cut applied to the SPI- 
DER sample (paper I). For Mb > — 19, ETGs exhibit a tight corre- 
lation between the metallicity gradient and either mass, luminosity, 
or log ctq. Brighter galaxies tend to have steeper gradients, as ex- 
pected by the lower efficiency of feedback processes in less bound 
(massive) systems (Larson 1974). At higher mass, colour gradi- 
ents exhibit a larger scatter, with no sharp depende nce on galaxy 
mass. It is currently not clear how the results of ISpolaor et al.l 
( l2009l) can be reconciled with the variation in the ratio of effec- 
tive radii with radius we find here. In fact, colour gradients are also 
determined by the change in the profile shape (i.e. the Sersic in- 
dex), besides radius, with waveband. Moreover, both metallicity 
and (small) age gradients can combin e to produce the observed 
internal colour gradients of ETGs (see iLa Barbera & de CarvaUiol 
1200^ ^. The tre nd of /r ^, , „ w ith r^^^p is consistent with a re- 
cent finding bv lRoche et al.l ( l2009h . who analyzed how the ratio of 
effective radii measured in g and r (using SDSS) correlate with sev- 
eral galaxy properties, for different families of ETGs (normal E/SO 
galaxies and BCGs). Although limited to the optical regime, they 
find that the mean ratio of radii measured in g and r become flatter 
for larger galaxies (Fig.|7). The trend of t'opt/'^n/j? 
plained by the increasing importance of dissipationless mergers in 
the formation of more massive galaxies with galaxy mass. Indeed, 
dry mergers are expected to wash ou t internal differences of stel- 
lar po pulation properties in galaxies l lWhitelll980l :l di Matteo et al.l 
I2OO9I) . A major role of dry mergers in the formation of mas- 
si ve ETGs has also been suggested, in a theoretical framewo rk, 
bv lNaab. Khochfar. & Burked feOOd) and lde Lucia et al.l ( I2OO6I) . 



10.3 The FP from g through K 

LBM08 derived the FP relation in the r (SDSS) and K (UKIDSS) 
wavebands, showing that the FP slopes exhibit only a small vari- 
ation with waveband, and that this variation is degenerate with 



respect to (i) the gradients of stellar population properties (i.e. 
age and metallicity) with galaxy mass, (5(log t)/5(log M) and 
5 (log Z)/S {log M), and (ii) the fraction of the FP tiU, /, which 
is caused by stellar populations. One main result of the present 
study is that using the grizYJHK coefficients of the FP we are 
able to break this degeneracy. The resulting probability distribu- 
tion of / is sharply peaked around zero, implying that the tilt of 
the FP in the NIR is not due to stellar populations. This result is 
in agreement with that of Trujillo et al. (2004), who found that the 
slope of the M/L vs. luminosity relation in K-ban d can be entirely 
due to structural non-homology of ETGs (see also lBusarello et all 
1 19971 : iGraham & Collesslll997l) . In B band, they found that a mi- 
nor, but still significant fraction (one-quarter) of the tilt is due to 
stellar populations. The re sults of Trujillo et al. (2004) contrast 
those o f iBolton et all ilOOH) . who argued that the tilt is more likely 
caused by a variation of the dark matter content with mass, with 
stellar populations pla ying a minor rol e, which fully agrees with 
our finding. Recently, I Jun & Inj j2008h have derived the FP rela- 
tion for a sample of fifty-six ETGs in the visible (V), NIR (K), and 
MIR (Spitzer IRAC) wavelengths and concluded that the slope "a" 
of the FP increases with the waveband. However, the uncertainties 
(see their tab. 2) seem to be still too large to conclude if "a" in- 
creases even further in the MIR wavebands. 

Spectroscopic studies of stellar population properties in 
ETGs have found that the (luminosity-weighted) age of ETGs 
tends to increase along the galaxy sequence, as parametrized in 
terms of either velocity dispersion or stellar a nd dynamical mass 
(e.g. iThomas et ai] |2005L iGallazzi et alj|2006h . The ages are usu- 
ally estimated comparing line spectral indices wit h the expecta- 
tions f rom stellar population models. In particular, IGallazzi et al.l 
l l2006i) found that the slope of the log t vs. log M relation is 
0.115 ± 0.056 (see their tab. 4). This value is estimated for a sam- 
ple of ETGs with a dynamical mass A/>lO^'^A/0, with a limit- 
ing magnitude comparable to that we adopt here in this work. The 
value of (5(logt)/(S(log A/) from Gallazzi et al. (2006) is signifi- 
cantly larger than what we obtain here, although still marginally 
consistent within 2-a (Tab. llOb . Moreover, we have to consider that 
age and metallicity values from spectroscopic studies always re- 
fer to the central galaxy region. Aperture corrections are based on 
measurements of line spectral indices for small samples of ETGs 
at 2 ~ and apply only to a relatively smal l radial range, with 
R < Re^ (see |j0rgenseij|l997l) . IGallazzi et al.l 12006) adopt a dif- 
ferent approach and instead of correcting the indices, test how the 
stellar population parameters vary with redshift, up to z ~ 0.12, 
for galaxies with similar physical properties (e.g. dynamical mass). 
The main drawback of this approach is that it relies on the assump- 
tion that spectral indices and their gradients do not evolve with red- 
shift. Considering the redshift range {z < 0.12), large galaxies 
are still observed only in a radial region of R<R^. The values of 
5{logt)/S{log M) and 5(log Z)/S{log M) we obtain from the FP 
analysis describe the total stellar population content of ETGs, as 
the photometric parameters entering the FP are defined in terms of 
the total galaxy luminosity of the 2D Sersic model. The information 
encoded in the FP is more similar to that provided by the colour- 
magnitude relation, where galaxy colours are usually measured 
within a larger aperture than that sampled by spectroscopic stud- 
ies. In fact, in agreement with our findings. iKodama et al. 
showed that the small redshift evolution of the CM relation implies 
that (i) all the (luminous) ETGs are equally old and (ii) more mas- 
sive galaxies are more metal rich than less massive systems. 

In the framew ork of the SAURON p roject, for a sample of 
twenty-five ETGs, ICappellari et~ZI ( I2OO6I) found that the varia- 
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tion of the dynamical M/L is well correlated with the Hp line- 
strength, implying that most of the tilt of the FP (i.e. the de- 
viation of FP coefficients from the Virial Theorem expectation 
under the assumption of homology and constant M/L) is in- 
deed due to galaxy age varying with mass. This resu lt apparently 
contra sts with findings of Trujillo et al. (2004) and iBolton et al.l 
( l2007h . and with our results, where both / and 5(logt)/5(logM) 
are consistent with zero. However, as als o noticed by LBM08, 
68% of the galaxies in the lCappellari et al] ( [2006) sample are fast 
rotat ors and 20% have low velocity disper sion (ff=60-85 km s 
~ ). Zaritskv. Gonzalez. & Zabludofd bOOeh and iD'Onofrio et al] 
have found that the FP of spheroidal systems depends 
on the covered range in m ass and velocity dispersion (see 
also I Graham & Guzman|[20o3 and references therein), with the tilt 
becoming larger (s maller "a") for galaxies in the low cro regime. 
IJeong et all ( |2009|) derived the NUV and FUV FP of thirty-four 
ETGs from the SAURON sample. They showed that the tilt is 
significantly affected by residual star formation in ETGs, mostly 
found at low ao (<100 km s~^) . Hence, the above men tioned dis- 
agreement with the findings of ICappellari et al" I ( l2006l) might be 
explained by the different range of velocity dispersion and differ- 
ent selection criteria of both samples. It is important to remember, 
as we have shown in Sec. 17. 31 that different subsamples of ETGs do 
not share the same FP relations. When binning the SPIDER sample 
according to Sersic index and axis ratio, we find that the tilt of the 
FP becomes larger (i.e. the slopes of the FP decrease) by a small 
but detectable amount for galaxies with higher n and larger b/a, 
with the effect being mainly due to a difference i n the "h" coeffi- 
cient o f the FP. The result for n is consistent with ID'Onofrio et al.l 
( l2008l) . who found that in the optical regime the "h" coefficient 
decreases significantly as the Sersic index increase s, while "a" is 
constant. However, one should notice that .D'Onofrio et al.l | |2008|) 
did not account for the fact that galaxies in different bins of n have 
different distributions in the space of the FP variables, and, as we 
show in Sec. 17.31 this might prevent a proper comparison of FP co- 
efficients. The fact that the variation of the FP tilt among galaxies 
with different n and b/a is similar from g through K suggests that 
it is more related to differences of galaxy properties (structural and 
dynamical), rather than to differences in the galaxy stellar popula- 
tion content. Kelson et al. (2000) derived the FP of 56 ellipticals, 
lenticulars, and early-type spirals in the cluster environment at red- 
shift z ~ 0.3. In agreement with JFK96, they found that the FPs 
of Es and SOs have consistent slopes. They also found that the FP 
of early-type spirals has a larger tilt (smaller "a") with respect to 
that of ETGs, likely due to a variation of the luminosity-weighted 
age with galaxy mass. This result might explain what we find when 
binning the SPIDER ETGs according to their optical-NIR colours 
and the discy/boxy parameter 04. Galaxies with bluer colours and 
more pronounced disc-like isophotes tend to have a more tilted FP 
(mainly because of a smaller "b"), with this effect smoothly disap- 
pearing from g through K. 
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Figure 21. Comparison of the variation of age (upper panel) and metallicity 
(lower panel) with stellar mass from the grizYJHK FP and the predictions of 
semi-analytical models of galaxy formation. Black circles and eiTor bars are 
the same as in fig. 6 of deL06, and represent median values of luminosity- 
weighted age and stellar metallicities of model elliptical galaxies. Black 
error bars link the upper and lower quartiles of the age and metallicity dis- 
tribution in a given bin of stellar mass. The inagnitude limit of the SPIDER 
sample corresponds to a stellar mass limit of ~ 2 X lO^O Mq, marked 
in the plot by the vertical dashed line. The blue circles are the peak val- 
ues of the distributions of luminosity-weighted ages and stellar inetallic- 
ities for model elliptical galaxies from the updated semi-analytical model 
of Wangetal. (2008), where the WMAP3 cosmology (rather WMAPl as 
in deL06) is adopted. The peak values are computed by the bi-weight es- 
timator Error bars denote 2ct uncertainties on peak values. Model ellipti- 
cals are selected as those objects in the semi-analytical model with a stellar 
mass fraction in the bulge larger than 80%, and colour index g — r > 0.5 
(consistent with the distribution of ETG's colours for the SPIDER sam- 
ple, see paper I). Ages refer to redshift 2 = 0, for both models. The red 
lines are the result of the analysis of Sec. |9] Their offset is arbitrarily cho- 
sen to match the models, while the slope's are obtained from the values 
of 5(log t) /(5(log M) and (5(log Z)/5(log M) reported in Tab.[To]for the 
BC03 model, with t = 12 Gyr and Z/Zq = 1 . 



10.4 Comparison to semi-analytical models of galaxy 
formation 

Explaining the stellar population properties of ETGs is a lingering 
problem for theories of galaxy formation and evolution. In the hier- 
archical scheme of galaxy formation, larger systems assemble their 
mass at later times. Hence, if star formation closely follows the 
mass assembling, one would naively expect more massive galax- 
ies to have younger stellar populations, in evident disagreement 
with (i) the red colours and old stellar populations characterizing 



the massive ETGs, and (ii) the o bserved bimodality of galaxies 
in the colour-rnagnitu de diagram dStrateva et aDl200lh . As shown 
bv iKauffmannI {m^ . semianalytical models (SAMs) of galaxy 
formation in the CDM framework can account for point (i) because 
more massive systems are indeed those forming stars at higher red- 
shift. For massive galaxies, the galaxy bimodality is also repro- 
duced, by preventing cooling flows in the centre of dark-matter 
haloes. This is achieved with som e ad-hoc recipe, like the feedback 
from AGNs dCroton et alj|200^ . As a result, both the luminosity- 
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weighted age and metallicit y of ETGs increase fr om lower to higher 
mass systems (see fig. 6 of Ide Lucia et aljfiood . hereafter deL06). 
Since the grizYJHK FP sets strong constraints on the variations 
of age and metallicity with galaxy mass, the natural question is if 
the amount of such variations can be accomodated in the frame- 
work of current models of galaxy formation. Fig. [21] compares 
the logarithmic variation of age and metallicity per decade in stel- 
lar mass, <5(log t)/S{log Ah) and <5(log Z)/S{\og that we in- 
fer from the FP □ (Sec. |9) and the expectation from SAMs. The 
plot shows the variation of the mean luminosity-weighted age and 
stellar metallicities as a function of galaxy stellar mass, i\f*, of 
model elliptical gal a xies for the SAM of deL06 (black circles), and 
that of IWang et al.l feOOS), where the latter model has been up- 
dated according to the WMAP3 cosmology. Interestingly, we see 
that current SAMs are actually able to match the results obtained 
from the analysis of the FP from g through K. Massive ETGs 
(> ~ 2 X 10^") have essentially coeval stellar populations, 
with more massive galaxies being slightly more metal rich, by a 
difference in metallicity of ~ 0.1 dex per decade in mass, than less 
massive systems. 

In a forthcoming paper, we will continue the analysis of scal- 
ing relations of ETGs, by presenting the dependence of the FP from 
g through K as a function of the environment where galaxies re- 
side, and discussing the implications for current models of galaxy 
formation and evolution. 



11 SUMMARY 

In this contribution, we present a thorough analysis of the FP of 
ETGs using a homogeneous dataset obtained in two wide-sky sur- 
veys (SDSS-DR7 and UKIDSS-LAS). As far as the FP derivation 
is concerned, we discuss fitting procedure, bias due to selection ef- 
fects, bias due to correlated errors on and < >e, and how to 
obtain meaningful FP coefficients. Below we summarise some of 
the main findings of this paper: 

1 - We examine the KR for all the wavebands available and find a 
smooth increase in slope from g (~ 3.44 ± 0.04) to K (~ 3.80 ± 
0.02), while the scatter seems to be independent of the waveband. 
Although the KR is just a projection of the FP relation, these results 
serve as a benchmark at the nearby Universe and will be essential 
for studies of ETGs at high redshift, for which not always large 
samples exist to probe the FP. In agreement the waveband variation 
of the KR slope, we find that the ratio of effective radii measured 
in g to that measured in K, ( ), decreases as re,K increases. 

2 - We measure the waveband dependence of the FP with unprece- 
dented accuracy. The trends of the FP coefficients, "a" and "b" 
(see Eq. [ij, with waveband are all very consistent regardless of 
the sample used (magnitude- or colour- selected). When using the 
log (Jo fit, we find that "a" is consistent, within 2-a, from r to 
K and for g band "a" differs significantly by ~3(t, while "b" is 
all very consistent. Using the orthogonal fit, however, "a" signifi- 
cantly varies by 12% from g through K and "6" does not change 
at all. 

3 - The analysis of the face-on and edge-on projections of the FP 
indicate, first of all, consistency with the results obtained when ex- 
amining the FJ and KR. Moreover, the scatter around the edge-on 

The quantities 5(log t)/(5(log M*) and S{\og Z)/5(log KL) are com- 
puted from the values of (5(log t)/<5(log M) and (5(log Z)/(S(log M) re- 
ported in Tab.[T0] and the relation (5(log M) = (5(log Af*)/(1 - 7x). that 
holds for / = (see Sec.|9). 



projection is about twice smaller than that of the face-on's, indicat- 
ing that the FP is more like a band rather than a plane. 

4 - We test the sensitivity of the FP solution to the veloc- 
ity dispersion measurement used, log ao (STARLIGHT) versus 
log cro(SDSS-DR6). Although these two measurements agree re- 
markably well, the value of "a" is systematically smaller when 
using the STARLIGHT values of log ctq, while "b" is insensitive to 
both measurements. Also, we find that the waveband dependence 
of the FP is the same regardless of the magnitude range used in the 
analysis. 

5 - The sample analysed is formed by ETGs covering a certain do- 
main in galaxy properties, like axis ratio (b/a), Sersic index (n), 
r ~ K colour, and 04. The FP slopes vary significantly for ETGs 
with different properties in the following way: ETGs with larger n 
have lower "b"; "a" is smaller in the NIR for the n > 6 subsam- 
ple, and in the optical both subsamples have similar "a"'s; The FP 
of round galaxies has smaller "a" (and smaller "b") than the FP 
obtained for lower &/a ETGs - the difference is more evident in the 
NIR. Also, boxy and bluer (r — K) ETGs exhibit an FP with lower 
"b", with this difference disappearing in the NIR wavebands. 

6 - Finally, we show that current Semi Analytical Models of galaxy 
formation match the results here obtained from the analysis of the 
FP tilt from g through K. This analysis implies that the NIR tilt 
of the FP is not due to stellar populations: massive ETGs have co- 
eval stellar populations, and are more metal rich than less massive 
systems. This is one of the crucial points of the FP study presented 
here. 



APPENDIX A: THE MLSO FIT 

We consider two random variables, X and Y, related by the linear 
model: 

Y^Pi+P2X, (Al) 

where pi and p2 are the offset and slope, respectively. We indicate 
as X and y the outputs of X and Y. Assuming that the y values are 
normally distributed along the orthogonal direction to the line, the 
probability of observing a given x and y pair is: 

P{r)dr = (27rao) • exp [-rV(2(7^)] dr, (A2) 

where r is the orthogonal residual, r = {y — pi — p2 ■ x) ■ 

( 1 + P2 ) and Oo IS the orthogonal scatter around the relation. 
In case where a selection cut is applied: 

2/ < ci + C2X, (A3) 
with ci and C2 assigned constants, Ea. lA2l modifies as follows: 

P{r)dr = K{pi,p2,ci,C2;x)-exp [-r'^/{2al)^ f{y-ci-C2x)dr,{A4) 

where the function / is equal to one when its argument 
is smaller than zero, and vanishes otherwise. The function 
K{p\,p2,c\,C2\x) is obtained by the normalization condition 
J P{r)dr = 1. If no selection cut is applied (/ = 1 identically), 

one obtains K = {2-Kal^ ^^'^ , and we recover Ea. lA2l In general, 
the K is given by: 

K = {2naiy'^^ • 2 ■ [1 + erf{t)]-^ , (A5) 

with t = [(ci -pi) + {c2-p2)x] /{V2aoy/l+pl), and erf 
denotes the error function. For a given sample of data-points, the 
likelyhood, L, can be written as 
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(A6) 



where both sums are performed over the entire dataset. In the case 
of the KR, one has y =< fi >c and x =log (Sec.|5]l. The mag- 
nitude cut can be written as < ^ >e < A/h™ + 38.56578 + 5 log re, 
where Af^m is the magnitude limit of the sample. This expression 
is identical to Eg. I A4I provided that ci = Mum + 38.56578 and 
C2 — 5. The MLSO coefficients of the KR are then obtained by 
minimizing the L with respect to pi, p2, and cto. 



APPENDIX B: MATCHING THE MAGNITUDE AND 
SURFACE BRIGHTNESS DISTRIBUTIONS OF ETG 
SAMPLES 

We consider the case where a set of n galaxy samples, with run- 
ning indices i = 1 to n, are given. In the case of Sec. 17.31 we 
have n = 2, and the two samples are obtained by splitting the 
magnitude-complete sample of ETGs according to a given galaxy 
parameter p. First, we select the sample with lowest sample size. 
For such sample, we define the minimum and maximum values of 
absolute magnitude, Mmin and Almax, and the minimum and max- 
imum values of <fj, >e, < >c,min and < jj, >c,max, respectively. 
We then construct a grid in the magnitude-< p >e plane, over 
the rectangular region from Mmin to Mmax, and < /x >e,min to 
< /X >e,m.ax- For 3 giveu cell k over the grid, we count the number 
of galaxies of each sample in that cell, rii^k - We take the minimum 
value of Wi.fc, n^, among all the given samples. For each sample, 
we then randomly extract galaxies whose magnitude and <fi>a 
values fall inside the given cell. This step is performed for all the 
cells in the grid. The procedure provides a subsample of galaxies 
from each input sample, with all subsamples having the same num- 
ber of galaxies and the same absolute magnitude and < fi>c dis- 
tributions. The mean values of Mmin and AI,nax, among the sub- 
samples analysed in Sec. 17. 31 amount to about —24.6 and —20.55, 
respectively, while the mean values of < /i >c,min and < /i >c,max 
amount to about 15.2 and 27.2 mag / arcsec? . The step sizes in M 
and < fj,>e are chosen to be 0.2 mag and 0.2 mag / arcsec? , re- 
spectively. This makes the number of galaxies in each cell of the 
grid to be smaller than 40. We verified that either reducing or in- 
creasing the bin size in a given direction by a factor of two does not 
change at all the results presented in Sec. 17. 31 
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